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ABSTRACT 

A  proposed  neutron-neutron  scattering  experiment  in  which  the  neutron 
source  would  be  colliding  beams  from  an  underground  nuclear  explosion  was 
investigated.   In  particular,  the  accuracy  to  be  expected  was  determined 
so  as  to  allow  evaluation  of  the  worth  of  the  experiment.   Phase  shift 
analysis  was  performed  for  Gaussian,  Exponential,  Yukawa,  and  Square 
Well  spherically  symmetric  potentials.  The  energy  range  considered  was 
20  kev  to  2.6  MeV.  Scattering  length  was  found  to  be  determined  to 
within  6%  and  effective  range  within  k0%   for  a  10%  scattering  cross  section 
accuracy.   Potential  depth  would  be  determined  to  at  least  .5%  and  range  to 
.3%. 
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1  .   Introduction 

A  low  energy  neutron-neutron  scattering  experiment  utilizing  a 
single  underground  nuclear  explosion  has  been  proposed  and  studied.  |_2,£] 
Preliminary  studies  indicate  the  scattering  cross  section  might  be  meas- 
ured to  an  accuracy  of  10%  from  20  kev  to  2MeV.  The  efforts  of  this 
paper  are  directed  to  a  numerical  investigation  of  the  expected  scatter- 
ing with  special  attention  given  to:   (a)  the  possible  existence  of  the 
dineutron;  (b)  the  effect  of  changes  in  effective  range  and  scattering 
length  upon  the  range  and  depth  of  the  spherically  symmetric  Gaussian, 
Exponential,  Yukawa,  and  Square  Well  Potentials;  (c)  the  effect  of  changes 
in  effective  range  and  scattering  length  upon  the  scattering  cross 
section  from  each  of  the  potential  shapes  in  (b) ;  and  (d)  a  comparison 
of  calculated  cross  sections  from  the  four  potential  shapes  in  (b)  with 
the  shape-independent  cross  section. 

Precise  determination  of  the  neutron-neutron  scattering  parameters 
would  be  important  for  a  number  of  reasons  as  stated  by  Moravcsik.  QZl 
These  reasons  include:   (1)  an  investigation  of  charge  Independence;  (2) 
an  investigation  of  charge  symmetry;  and  (3)  the  existence  or  non-existence 
of  a  slightly  unbound  dineutron  resonance  state. 

This  investigation  attempted  to  locate  a  dineutron  state  in  section 

2  to  provide  a  theoretical  justification  for  its  experimental  quest.   In 
Sections  3  and  k,    2   to  16%  deviations  in  the  ratio  of  effective  range  to 
scattering  length  from  the  accepted  n-p  ratio  were  investigated  for  their 
effect  on  scattering  cross  sections  and  Gaussian,  Exponential,  Yukawa, 
and  Square  Well  Potential  parameters.  These  results  should  allow  some 
determination  of  the  accuracy  in  experimental  results  required  to 
determine  to  a  given  accuracy  the  shape-independent  scattering  parameters. 
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In  Section  5,  the  cross  section  from  effective-range  theory  is  compared 
with  cross  sections  from  potential  scattering  to  test  the  validity  of 
the  shape  independent  approximation.  This  report  concludes  with  a 
summary  in  section  6  and  the  complete  computer  program  used  in  the 
calculations  along  with  the  calculation  results  as  a  part  of  the  appendix. 

The  analysis  was  performed  using  phase  shift  scattering  methods. 
The  well  established  scattering  length  and  effective  range  for  neutron- 
proton  scattering,  a*  and  rQ  s*,  were  taken  as  a  base  for  calculations. 
The  ratio  of  effective  range  to  scattering  length  was  formed  to  yield 
a  spherically  symmetric  potential  range  and  depth  for  square  well, 
exponential,  gaussian,  and  yukawa  potentials.  The  absolute  value  of 
the  scattering  length  was  then  varied  in  2%   steps  up  to  16%  while 
holding  the  effective  range  constant  at  its  accepted  value.   17  sets 
of  potential  parameters  were  thus  obtained  for  each  of  the  four  poten- 
tial shapes.  Then,  using  the  maximum  and  minimum  acceptable  value  of 
effective  range  as  a  base,  the  effective  range  was  varied  in  2%   incre- 
ments up  to  6%  while  holding  the  scattering  length  constant  at  its 
accepted  value.  Eight  additional  sets  of  potential  parameters  were 
obtained  for  each  of  the  potential  shapes. 

The  spherically  symmetric  radial  wave  equation  was  then  numerically 
integrated  using  a  fourth  order  runge  kutta  method  for  all  calculated 
sets  of  parameters  and  potential  shapes.   Integration  was  performed  for 
30  values  of  center  of  mass  energy  ranging  from  20  kev  to  2.6MeV  and 
was  carried  out  to  a  radius  at  which  the  effect  of  the  nuclear  potential 
was  no  longer  felt. 

The  logrithmic  derivative  of  the  numerically  integrated  solution 
was  then  compared  to  the  logrithmic  derivative  of  the  exact  field-free 
asymptotic  solution  outside  the  range  of  the  nuclear  potential  and  a 


phase  shift  obtained. 

The  total  cross  section  as  a  function  of  center  of  mass  energy  was 
then  calculated  from  the  phase  shift  for  each  case.  A  ratio  of  each 
cross  section  to  that  obtained  from  the  accepted  neutron-proton  scatt- 
ering potential  parameters  was  formed.  This  gave  the  per  cent  change 
In  total  cross  section  for  each  incremental  change  in  scattering  length 
and  effective  range. 

Per  cent  change  in  total  cross  section  from  that  for  n-p  scattering 
parameter  calculations  was  plotted  against  energy  for  significant  incre- 
ments in  two  different  energy  ranges,  20  kev  to  100  kev  and  100  kev  to 
2.6MeV.   Potential  depth  and  range  were  plotted  against  the  ratio  of 
effective  range  to  scattering  length  for  each  potential  shape.  A  com- 
parison of  shape- independent  effective  range  scattering  calculations 
with  potential  scattering  calculations  was  made  by  plotting  the  total 
cross  section  for  each  case  on  the  same  graph. 


2.   Possible  dineutron  state 

An  attempt  was  made  to  justify  the  possible  existence  of  a  dineutron 
resonance  state  near  zero  energy  using  a  compound  nucleus  model,  but 
failure  was  met  in  attempting  to  formulate  a  compound  nucleus  model  for 
neutron-neutron  scattering. 

Numerical  calculations  for  neutron-neutron  potential  scattering 
were  made  for  a  large  number  of  potential  depths  on  the  Control  Data 
Corporation  1604  digital  computer  in  an  attempt  to  find  possible 
resonances,  particularly  a  sharp  resonance.  Gaussian,  exponential, 
and  yukawa  potentials  were  used  with  potential  depths  ranging  from 
0  to  the  value  obtained  from  the  accepted  n-p  scattering  parameters 
in  2.5  MeV  steps  over  two  center  of  mass  energy  ranges,  20  to  100  kev 
and  .1  to  2.6  MeV.  Results  were  scanned  for  resonances.  No  reson- 
ances were  expected  and  none  were  found. 


3.  Effective  range  and  potential  parameters 

Detailed  comparisons  with  experiments  have  given  for  the  singlet 
effective  range  and  scattering  length  for  neutron-proton  scattering: 

r0*^  2.+  6±  0.12   FM  a*s  -  -  23  .  677   ±  0.02<*  FM      L4CI 

For  an  examination  of  the  effect  of  deviations  from  this  scattering 
length  on  the  scattering  cross  section,  the  ratio  ra>s  /  <x   was  formed 
and  a  was  varied  in  magnitude  from  0-16%  in  2%  steps  while  holding  r0)S 
at  its  accepted  value,  2.46  Fm.  To  examine  the  effect  of  deviations  in 
effective  range  on  the  scattering  cross  section,  the  ratio  r0)%    I  clx 
was  again  formed  using  the  accepted  value  for  as,  -23.677  Fm,  and  the 
maximum  and  minimum  values  for  roS      ,2.58  Fm  and  2.3*+  Fm  respectively. 
ro s  was  varied  UP  to  6%  in  2%   steps  in  the  direction  of  increasing 
deviation  from  the  accepted  value,  2.46  Fm. 

Four  spherically  symmetric  potential  shapes  were  then  examined. 
These  potentials  are: 

Square  Well:    V  •  "  V0     (~  £  rn 

v  -  o  r  >  rn 

Gaussian  Wei  1 : 

V  =  -Va  ^K-t/rn)*] 

Yukawa  Well: 


Exponential  Well:  "• 


v  =  -  Vo  e*f  L-  r/ml 


Blatt  and  Jackson  have  arrived  at  a  correspondence  between  the 
ratio  of  effective  range  to  scattering  length  and  the  potential  para- 
meters, V0  and  r  ,  for  each  of  the  well  shapes  investigated,  [j] 
The  well  parameters  are  given  in  terms  of  the  intrinsic  range,  b,  and 
the  well  depth  parameter,  S.  The  well  depth,  V/r\,  for  the  potential 
in  question  is  multiplied  by  a  constant  such  that  the  first  resonance 
occurs  at  zero  energy.  The  well  depth  parameter  Is  then  the  inverse 
of  this  constant  given  by: 

V(r>=  5V„s(r) 

The  intrinsic  range,  b,  is  determined  by  the  shape  and  range  of  the 
adjusted  potential  through  the  radial  wave  function  for  that  potential, 
It  is  defined  by: 

RES       o  >6S 

Blatt  and  Jackson  provide  interpolation  formulas  to  relate  VQ  and 

rn  to  rQ  and  a.  The  form  of  their  formulas  have  been  changed  slightly 

to  fit  the  definition  of  the  potentials  in  this  report.   In  their 

modified  form,  these  interpolation  formulas  are: 

Square  Wei  1 : 

S  =  Ll-  I.  SWr./o)* .353(r%/fc}a]/[j-  E.3M-  (r./o.)  t  1.3  0  (To/oO^] 

\>  =  W C I  -.?O0(s.-|)-h4-86  (s- D*-.4 76 (*-!)*  +  .376(s- 1)4"-.  |  5"6  Cs-l>S  J 

V0=  1.01176  XI  0     S/fc£ 

rn-\> 

Gaussian  Wei  1 : 

s*  Ll-2.030  UVa.)  4.5e7  (r0/(L)1  3 / Li  -  a.eio(r6/0L.)4l.fc0ff  cr0/o>z3 

fc»  =  r6/ci-.fe4us-o  +.5"6fc<s-i)a'-.7£i  c^-i^-vi.nfccs-o^-i.ois-fs-if  J 

Vo*  c.aqeos  x  io2"   s/b2 


Exponential   Well: 

s  =  ci-fc.<m<Ur0/o^  +  i.3  6oa(re/o.t3  /U-  3.rm4>{r0/oO  +  'b.m7Cn,/<OJ 

b*  rft/[.|-.qo4C?.-0  +.145Cs-l^-.5-43C*->'>  +  .H4  C^-lV^ 

Vo-  7.  SI  5-4  I  XIO      s/b* 

rn  z  t/3.r4^ 

Yukawa  Wei  1 : 

S^Cl-  3.&oS3<ro/0L>+Z.0  SS-^ro/o^J/Cl-  4.  \  34fc  <.ro/0w>  + 'i.ioiaoWo&J 

V»b  3.1*9 SLUM  KI6*14&/W*- 

The  results  of  these  calculations  are  presented  both  In  tabular  and 
In  graphical  form  In  tables  I  and  II  and  In  graphs  1  through  16.  The 
tables  are  arranged  to  exhibit  V0  and  rn  as  a  function  of  the  percentage 
deviation  of  r  and  a  from  r*    and  a*s,  the  best  values  from  n-p 
scattering.  The  graphs  depict  V0  and  rn  as  a  function  of  r0/a  using 
the  left  and  bottom  scales  respectively.  Using  the  top  and  right 
scales,  the  percentage  deviation  of  VQ  and  rn  from  the  best  n-p  scatter- 
ing parameters  Is  shown  as  a  function  of  the  percentage  deviation  of 
r0/a  from  r*0fS/a*s> 

Variations  In  potential  parameters  were  found  to  be  of  the  order 
of  a  .5%  Increase  In   |V0|  and  a  .5%  decrease  In  rn  for  a  10%  increase 
in  scattering  length.  A  10%  increase  in  effective  range,  however,  was 
observed  to  cause  a  decrease  in   I  Vol  of  approximately  8%  while  increasing 
rn  about  9%.  This  data  is  shown  in  table  III. 


< 
3 


vOJ1  M-  O00  NIA4-  N  -  CT\  00  vo  LO  J"  CM 
0O0O0O0O00  rvfv.(Nh»NNvOvOvO\0\0\D 

ooooooooooooooooo 
c        

I.  —    —    —   —    —    —    —   —   —    —    —    —    —    —    —   —    — 

cmc-vj-iavo  r»»oo  o^  o  o  —  n  roj-  iavo  r» 
o       • 

I  UM/MAUMAUMAUMnUMAlAUMAlAUMn 


to 

LU 

I- 

< 


< 

i 


o 
o. 
x 

uJ 


c 
u 


inoO  Or*Mf\CO  O  CMAOO  -  (AvO  ff\(M  tfiCO 
vO  \D  vO  vD  vO  vO  vO  vO  vO  vO  vO  vO  vO  \0  vO  vO  vD 


NwinvONcno-Nj-  u\\o  oo  <r»  o  -  o 

O       OOOOOO  —  ——  —  —  —  —  —  —  CMCMCM 

>    CO  CO  CO  CO  CO  CO  CO  CO  CO  CO  CO  CO  CO  CO  CO  CO  CO 

I         ————————————————— 


UJ 

3 

c 

1- 

l_ 

o 

z 

a. 

< 

z 

CO 

o 

CO 

a. 

=> 

3 

< 

o 

p— 

o 

> 

X 

1 

LU 

H 

_J 

o 

00 

z 

< 

UJ 

U\J-  CM  —  CT\00nO  LO  CO  CM  O  CT>00nX>  ltvcocm 
-$-^-$-£cocococococococmcmcmcmcmcm 
vo  *«D  vD  vD  vO  vD  vO  vO  vO  vO  vO  vO  vO  vD  vO  vO  vO 


00<n-N4vOrNO\ON4-lANCO   O  CM  CO 

u\lo\OvO\DvOvOvo  f**.r«»r*-.r,~>r,^r>.oooooo 
cooooooocooooooooooococooooooooooo 

CO  CO  CO  CO  CO  CO  CO  CO  CO  CO  CO  CO  CO  CO  CO  cO  CO 


o 


< 

o 

CO 


CO 
LU 

o 

z 
< 

X 

o 


o 

LU 

u. 


LU 

3 


a: 

I 

CO 


CTvr^LOCOCs!  O  CO  VO  LOCO—  CTVOO  vo  J"  ^M  — 

rvNrs.h.NrsvovovDvOvO  lololololola 

CO  CO  CO  CO  CO  CO  CO  CO  CO  CO  CO  CO  CO  CO  CO  CO  CO 

cmcmcmcmcmcmcmcmcmcmcmcmcmcmcmcmcm 


CMCViCMCMCMCMCMCMCMCMCMCOCOCOCOCOCO 
0\  Q\  CT\  CT\  O.  CT\  C7\  CT>  <3\  0\  OS  OS  OS  OS  OS  OS  OS 

v0v0vovOvovOvov0vov0vovov0v0vov0vO 


vO 

CM 
II 

O 

L. 


LOCO  —  COvO-tf'OMOOOvO-TCMOOOvo.d'CM 

Mnmooovfij  cM<nr«»tnco  —  oovO-d'  cm 

<D  CMCO-3-LOLOVDf**OOOO0\O-   N  N  «J"   U\ 

\       r».cT\— coLor-tTv«— colocoo«mJ"vOooo 

O       0000vT»CT\CTVCT\CT*OOOO  —  —  —  —  —  CM 
«-         OOOOOOO  —  —  —  —  —  —  —  —  —  — 


I     »     I 


I 


( 


« 


•     I     I 


I     I 


4) 

* 

ID 

<0 

C 

(D 

c 

JO 

<_ 

& 

vD  J-  NOOOXDJ1  N         N  J-vOOO  O  NJ-vO 
I      I      I      I      I      I      I      I    O 


2 


o 
<*\^  —      .4"  ca  cm  .^  la  la  ca      —  cn  cn 

U\vi  N      VO  j-  nnnooo        r-oioo 

VOJ-N          —  —  —  O   O  O  ON         CO  vO  J" 
C vOJ- 

L.  _   —   ^_  _   __    _   _    _   _  .. 


-d-  rA-d"  00  OOO  m-  U\  CO  — 

N  O   N  N-   OO   rsNOO  LA      »      • 

0      ■  •    •     ♦  -  ••     ••<•     •     •  •  c*nI  r*^ 

>       noo  oo  oNNNtns  smm 

l         Nmj  «a-  la  la  la  vO  vO  vO  ca  .—  ca 


—I  —  Cn  -3*  CO  —  CA  LA  la  ^  vO  00  vO  —  00  vO 

LU  p-m4M0OlANN  NOOOOOO  OO 

5  o<m^«-  o  cnoo  la  cm  .—  0>CAN  O  N 

</>  C  •COI^Nf^vOvOvOvOvOlAUMAJ'N 

or  -j       »-       — 

Ul  < 

s       i 

<  iu  -j-Nooo-vONr^inNNo-oo 

O?  SC  NNO      •    ' N     • 

<  o         o       •    ••cnoLAcn.—  j'ovo  —  cA»-d* 

CL  O-  >  N00  NO-  —  —■  W  J1  UMAVO  O  lAO^ 

X  I  LAvO  CA  ■—  —  —  —  —  •—  —  —  —  NMN 

-I  UJ 

< 

i 

Ul 

o 

a.  j  r>.  r*.  ca  cn  r»»  on  cm  cn 

r»  — co  — 

O  2K  _C-      •_"■_•'*     *  *      •  CAvO 

Ow 


-J 

r>».  r*.  ca  cn 

-J 

j--a-  -a- oo 

Ul 

IAN  c\rs 

■» 

c 

.       .       •       • 

l_ 

CM  CM  —  •— 

z 

< 

MM 

(/) 

OO 

o  —  en  la 

3 

ca  en  on  o 

< 

O 

.... 

o 

> 

LA  CAvO  CM 

1 

—  —  CM  CA 

oo  cm  vo  CO 
<  o  ca  cn  on  o  vocm 

^         0£  ^  O  •     •     *     •  •     •  vO  00 

r»*  o  o  c% 

■4"  vo  —  CM 
uj       > 

_i       — 

CO  t- 

<         O 

u. 

Ul 

CA  CA  r«-         —  CA  LAIA  CA  O  VO  vO   —   O 

z  _i  o  vo  —       r».  cm  r^vo  la  —  vo       o  -d-  oo 

—  _l  1*^  CM  CO  LA  LA  J"  CA  CM  CM  —  Cn  J"  vO 

1 1 1         c       ...        .......        ..  en 

L.  CACACM  CMCMCMCMCMCMCM  —  —      • 


Nj'vo  rAcncncM       r-»  en  nrxrv 

N4  n  cm  r*» ca en  r»»  j-  cm  j- en    • 

O       vO  vO      *  .......  ..  LA 

>        •    .  —  ^r  -a-  lavo  oo  en  o  vo  vo  o 

I           VOOO—  —  —  —  —   —  —  CM  CM.4-  — 


r^oooooo  lavo  r^cAocAr^cncncn 
calai^oo  cMj-vocncAiAr*-o  —  cao 
cm  j-  vo  cm  ca  —  cnoo  oo  oo  oo  LA  —  CA  vo 

VOLA^J-CA  —  00CAC0vO-*CACACMLA 

id      voj-cm  —  —  —  oocncnoNcnoovo  — 

^ . -  —  —  —OOOOOOJ" 

o o 

1.        I    I     I    I     I     I     I    I     I     I     I     I     I     I      . 


u> 

5 

Ul 

o 

Ul 

z 

ec 

S 

< 
=3 

o 

r»» 

O* 

r>» 

</) 

u. 

vO 

o 

• 

fA 

t- 

<M 

o 

1 

Ul 

u. 

m 

u. 

M 

Ul 

<D 

V 

cn   o 

C   «- 

00  CO  00         00  00  CO 

10 

O   O   O   O  00  0O  CO         CO  CO  00 

-C    c 

vO  -4"  CM  —     •     •     «            «      •      'OOOO 

8- 

CO  VO  J"         J- vO  CO  —  <M  J- vO 

O     1      1      1      1      1      1      1 

TABLE  III 
PERCENTAGE  CHANGE  IN  WELL  PARAMETERS  CAUSED  BY  A  TEN  PERCENT  INCREASE 
IN  SCATTERING  LENGTH  AND  EFFECTIVE  RANGE 
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k.      Phase  shift  analysis 

A  partial  wave  analysis  of  neutron-neutron  scattering  was  performed 
for  all  potential  parameters  calculated  in  section  3. 

Assuming  a  spherically  symmetric  potential,  the  radial  wave  equa- 
tion 

*»!r(rt£)  +  (^i>-V«->l  -  ^]  R=  O  (1) 

becomes 


-~k ajH  +o<^  *  m^n  u  =  eul  (2) 

where    H(r)--rR(r) 
Rearranging,  equation  (2)  becomes 

£Ut  *-{**■- v'crt-  41^3  U  (3) 

where      **«  NE/ti*   -  [E(MEV>/4I.4<KI  *lOU 

and     V(r)  aMV(rt/hfc»cycrt(KfeO/4l.4^XlOat         Q6] 

Ms  neutron-proton  reduced  mass 
E=  center  of  mass  energy  =  1/2  laboratory  energy 
Separating  equation  (3)  into  two  first  order  equations  for  simul 
taneous  numerical  integration  gives 


ft— ^-v'cr>-4itop]U 


(5) 
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Equations  (4)  and  (5)  were  numerically  integrated  using  a  fourth 
order  Runge-Kutta  method  with  Boundry  conditions: 
For  r  =  O 

UCo)  -O 

Z  (CO  =  ± 

Integration  was  carried  outside  the  range  of  the  nuclear  potentials 
with  500  to  750  steps  to  about  15  fm  for  two  center  of  mass  energy 
ranges;  0-100  kev  in  20  kev  intervals,  and  .1-2.6  HeV  in  .1  MeV  intervals 
The  values  of  Z  and  U  at  r  -  15  fm  are  referred  to  as  Zn  and  Un. 

The  exact  asymptotic  solution  of  the  field-free  wave  equation, 
equation  (1)  with  V1  -  0,  is  given  by  Zeleny  as:  Q7J 

RjCO^Asj  ^cos  gA  J  ikr}  -sin  &,  n;  (fcrtj  (6) 

from  which 

With  the  substitutions 

equations  (6)  and  (7)  become: 

lyo*  Ajr{eos  *,  JjCkrt  -  &in    s^tfe^}  (8) 
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where  6.  fs  the  scattering  phase  shift  and  j«   and  n   are  spherical 

Bessel  and  Neumman  functions  respectively. 

The  scattering  phase  shift,  g$    ,  i s  obtained  by  matching  the 

logrithmic  deviations  of  the  exact  field-free  and  numerical  solutions 

at  a  range  where  the  nuclear  force  is  negligible.  This  procedure 

yields ; 

For  X   =  0 
Znl   ^ 


=  i.rLM  n,(kr)  sin  &0  -  k\x(kr)  cos  &o  +  j0(fer>  cos  s< 


nispw 


■-n.(krt  sin  4»3/tfti<feit  cos  &*  -  n0tfcrt  sm  s^     \r^ 
For    y  z  1 


rsisFw 


+  cosS1i,(fer>-s»nSln,(*?r5J/£rtos&,j((»?rV^n  S^Cter^j 


5  FM 


Rearranging  and  taking  the  inverse  tangent  we  get  the  following  express- 
ions for  the  phase  shift: 

§.  --  HhaLn'ffifc  4ol»rVr-*0,Cter>]r  -  a.  Up") 

i(^Anc<»ir^  +  ^n,cJ*e)]r  -  no(fer)  J 


(10) 


§,  - 1  an 


(ID 


The  total  cross  section  as  a  function  of  energy  is  obtained  from 
the  phase  shift  using  the  relation  [$\ 

<?CO    =    Z  if[   Ce.v+0  ^n?^ 


(12) 
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For  Identical  particle  scattering,  equation  ( 1 1)   becomes,  for  the 
singlet  state: 


even 
and  for  the  triplet  state: 


«i(E)*  £  irr  Ce.tf+0  fc»n*<Sjj 


% 


The  computation  of  cross  sections  was  made  for  both  the  singlet 
1=0  state  and  the  triplet  X   =  '  state.  The  j(  =  1  cross  sections 
were  observed  to  be  at  least  3  orders  of  magnitude  smaller  In  the 
energy  range  under  consideration  and  were  disregarded  in  subsequent 
computations.  The  computer  program  in  appendix  III  however,  retains 
the  capability  of  X   ■  1  calculations. 

The  computer  program  In  appendix  III  was  used  in  cross  section 
calculations  for  all  the  potential  parameters  in  section  3.  First 
the  program  computed  the  potential  parameters  at  the  base  scattering 
length  and  effective  range  and  at  percentage  deviations  from  this  base 
read  Into  it  in  its  input  section.  The  total  scattering  cross  section 
was  then  computed  as  a  function  of  center  of  mass  energy.  The  calcul- 
ated cross  section  for  each  set  of  potential  parameters  was  then  divided 
by  the  cross  section  for  the  base  parameters  to  normalize  the  cross 
section  and  present  it  as  a  deviation  from  the  base. 

The  resulting  data  was  plotted  to  show  the  percentage  difference 
in  total  cross  section  from  the  calculated  total  cross  section  for  n-p 
scattering  with  ag  ■  -  23.677  fm  and  r0>s  s  2.46  fm  as  a  function  of 
energy  for  percentage  deviations  of  a  and  r  .  The  graphs  are  drawn 
over  two  distinct  energy  ranges,  20  -  100  kev  and  .1-2.6  MeV,  to 
better  portray  results  and  separate  graphs  are  provided  for  each 
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potential  shape. 

In  Moravcsfk's  paper,  the  function  0_(T)  Is  defined  by: 

«">  -(H)"' a  k+  M "  £)  fr  s  4TT/«^ 

where  e%MT  for  T  ■  laboratory  kinetic  energy  and  m  z   neutron  mass. 

2. 
Moravcsik  has  determined  that  in  a  10%  experiment,  a  linear  approxi- 
mation for  0_(T)  would  be  sufficient. 

The  value  of  this  straight  line  at  £  s  0  would  give  l/a2  directly 
while  the  slope  of  the  straight  line  would  determine  rQ/  and  the  sign 
of  a. 

Moravcsik  concludes  that  a  determination  of  total  cross  section 
to  10%  accuracy  at  ten  suitably  distributed  energy  points  could  det- 
ermine the  scattering  length  to  3%  and  the  effective  range  to  50  -  70%. 

Tables  IV,  V  of  this  report  show  that  a  10%  determination  of 
total  cross  section  in  the  vicinity  of  the  n-p  scattering  parameters 
will  allow  a  6%  determination  of  as  at  20  kev  and  a  16%  determination 
at  400  kev  while  at  2.6  MeV  a  16%  error  in  scattering  length  has  only 
a  2%  effect  on  total  cross  section.  The  same  tables  show  that  a  10% 
determination  of  total  cross  section  will  allow  a  k0%   determination 
in  effective  range  at  2.6MeV  while  a  100%  error  in  effective  range 
has  approximately  a  2%  effect  at  20  kev. 

An  alternate  method  of  determining  the  effective  range  parameters 
is  therefore  offered.  The  magnitude  of  the  total  cross  section  at 
20  kev  could  be  compared  to  supplemented  data  of  this  report  to  arrive 
at  a  value  for  a  within  6%.  Graph  17  is  shown  as  an  illustration  of 
possible  supplemental  data.  The  magnitude  of  the  total  cross  section 
at  2MeV  or  greater  could,  in  the  same  manner,  be  used  to  arrive  at  a 


31 


TABLE   IV 
EFFECT  OF  A  10%  ERROR  IN  TOTAL  CROSS  SECTION  UPON  THE  ERROR   IN  a. 


%   Ei 

rror  In  a 

for 

ENERGY (kev) 

SQUARE  WELL 

GAUSSIAN 

EXPONENTIAL 

YUKAWA 

20 

8-16 

6 

6 

6.5 

60 

8-16 

8.5 

8.5 

9.5 

100 

8-16 

11 

12 

16 

200 

8-16 

16.5 

16.5 

18 
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TABLE  V 
EFFECT  OF  A  10%  ERROR  IN  TOTAL  CROSS  SECTIONS  UPON  THE  ERROR  IN  r( 

%error  In  r_   for 


ENERGY 

(MeV) 

GAUSSIAN 

EXPONENTIAL 

YUI 

1.6 

60 

58 

62 

1.8 

56 

52 

59 

2.0 

53 

50 

57 

2.2 

50 

k7 

55 

2A 

hS 

k5 

52 

2.6 

ko 

hi 

50 
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GRAPH   17 


value  for  rn   c  within  50%.  Because  of  the  negligible  effect  of  both 
0,5 

large  errors  In  rQ  s  at  the  low  energy  cross  section  and  large  errors 
in  as  at  the  high  energy  cross  section,  they  may  be  determined  indepen- 
dently, each  having  little  effect  upon  the  determination  of  the  other. 
Either  gauss i an,  exponential,  or  yukawa  data  could  be  used  because  of 
only  about  a  1%  difference  in  as  value  among  them  for  any  given  cross 
section.  Thus  the  scattering  length  and  effective  range  could  be  det- 
ermined with  a  minimum  of  two  points,  one  at  either  end  of  the  energy 
range.   For  greater  accuracy,  the  computer  program  in  the  appendix 
could  easily  be  adapted  to  a  searching  interpolation.  An  advantage 
could  be  gained  if  the  necessity  for  only  two  values  of  cross  section 
would  allow  a  more  precise  determination  of  those  two  values.   It  can 
be  seen  from  accompanying  graphs  that  the  energy  must  be  accurate  to 
10%  to  achieve  meaningful  results  In  this  way.   If  this  energy  accuracy 
Is  not  feasible,  then  the  method  outlined  by  Moravcslk  for  scattering 
length  and  effective  range  determination  must  be  used.   In  any  event, 
Moravcsik's  development  of  expected  accuracies  is  confirmed  by  this 
phase  shift  analysis. 
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5.  Shape  independent  approximation 

The  total  cross  section  for  neutron-neutron  scattering  was  cal- 
culated by  the  shape  independent  approximation  of  effective  range 
scattering  theory  using  r*Q  s  and  a*  .  This  cross  section  was  compared 
with  that  from  the  potential  scattering  calculations  in  section  k. 

Segre  gives  for  the  total  cross  section 

where  a(k)  is  the  scattering  length  as  a  function  of  wave  number.  An 
approximation  for  a(k)  is  given  for  k  ^  0  as: 

which,  when  substituted  for  a(k)  in  the  cross  section  formula  yields: 

Formula  (la)  was  used  in  a  calculation  of  total  cross  section  in 
the  energy  range  .1  -  2.6MeV  and  compared  with  the  cross  sections 
obtained  from  potential  scattering  for  the  same  parameters.  The  results 
are  shown  in  graphs  32  -  33. 

The  results  of  the  effective  range  theory  calculations  compare 
favorably  with  potential  scattering  calculations.  A  2%   difference  is 
the  greatest  deviation  noted  in  the  comparison.  This  confirms  the 
commonly  held  belief  that  the  shape  independent  approximation  of 
effective  range  theory  is  quite  accurate  in  this  energy  range. 
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6.  Summary 

In  section  k   it  was  determined  that  a  10%  scattering  experiment 
could  define  a  scattering  length  to  approximately  6%  and  an  effective 
range  to  approximately  40%  accuracy.  This  accuracy  might  be  achieved 
with  only  two  values  of  cross  section;  that  at  20  kev  and  that  at  2.6 
MeV,  if  the  energy  could  be  determined  to  10%  in  both  regions.  Certainly 
It  would  be  achieved  by  a  least  squares  fit  at  10  -  20  appropriately 
distributed  energies  as  proposed  by  Horavcsik. 

In  the  region  of  the  accepted  n-p  scattering  length  and  effective 
range,  it  was  found  that  a  2%  change  in  scattering  length  resulted  In 
a  3.1%  change  in  cross  section  at  20  kev  center  of  mass  energy  while 
only  causing  a  .1%  change  in  total  cross  section  at  2.6MeV.  A  20% 
change  in  effective  range  was  seen  to  cause  a  5%  change  in  cross  section 
at  2.6MeV  while  having  almost  no  effect  at  20  kev.  Thus  the  conclusion 
was  reached  that  effective  range  and  scattering  length  could  be  determined 
separately  and  independently  by  the  cross  sections  at  the  high  and  low 
ends  of  the  proposed  energy  range. 

In  section  3  it  was  found  that  a  1%  change  in  potential  depth 
or  range  resulted  in  about  a  20%  change  in  scattering  length  and  about 
a  1%  change  in  effective  range  for  Exponential  and  Gaussian  wells. 

For  the  Yukawa  well,  a  1%  change  In  potential  depth  and  range 
resulted  in  about  a  12%  change  in  scattering  length  and  a  1%  change 
In  effective  range. 

It  was  conversely  determined  that  a  6%  determination  of  scattering 
length  can  determine  the  Gaussian  and  Exponential  potential  depths  to 
less  than  -   .2%  and  the  potential  range  to  less  than  -  .3%  for  the 
n-p  scattering  length  and  effective  range  vicinity.  The  Yukawa  well 
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depth  would  be  determined  to  less  than  -.5%  and  the  potential  rat 
to  -.2%. 

It  therefore  seems  quite  likely  that  a  10%  scattering  experiment 
would  determine  the  scattering  length  and  effective  range  accurately 
enough  to  furnish  an  acceptable  set  of  potential  parameters.  These 
potential  parameters  appear  to  be  adequate  to  investigate  charge  - 
independent  and  charge-symmetry  effects. 

Effective  range  theory  was  quite  adequate  to  determine  cross 
sections  in  the  examined  energy  range.  A  1%  agreement  was  found  with 
potential  scattering.  There  should  be  little  question  of  i£s  propriety 
here. 
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APPENDIX  1 
COMPUTER  PROGRAM  DESCRIPTION 


The  program  Is  written  in  FORTRAN  II  language  and,  as  listed  In 
Appendix  III,  is  arranged  to  operate  with  the  Control  Data  Corporation- 
1604  as  a  central  processor  and  the  IBM-  1401  as  peripheral!  equipment. 
The  library  tape  functions  used  are  the  standard  FORTRAN  0  1  functions 
with  the  exception  of  a  subroutine  called  DRAW.  This  subroutine  is  a 
part  of  the  library  tape  at  the  Computer  Center  of  the  United  States 
Naval  Postgraduate  School  and  is  identified  by  its  large  computer 
operators  co-op  ID:  J7=NPS-DRAW.  The  subroutine  was  written  by  J.R. 
Ward  of  the  U.S.  Naval  Postgraduate  School.  Where  this  subroutine  is 
not  available,  the  program  can  be  run  by  omitting  subroutines  SKETCH 
and  NOW  and  their  call  cards. 

The  program  consists  of  a  main  routine  and  13  subroutines.  The 
main  routine  calls  the  subroutines  in  proper  sequence  and  also  provides 
for  a  sequential  variation  of  the  parameters  rn,  and  V    .     The  subroutines 
and  the  function  each  performs  are  listed  as  follows: 
Subroutine  INPUT2      Reads  in  data  cards 

Calculates  and  tabulates  energies,  wave  numbers, 
scattering  lengths,  effective  ranges,  potential 
parameters,  and  calls  bessel  functions. 
Prepares  for  and  calls  integration  routine,  forms 
the  logrithmic  derivative  of  the  integrated  function, 
Integrates  equation  provided  by  DER5V  for  f  s  0, 
Integrates  equation  provided  by  DERflVl  f©rj(s  1. 
Provides  J(s  0  equation  for  integration  for  a 
particular  potential. 
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Subroutine  READY2 


Subroutine  FIGURE 

Subroutine  RKUTTA 
Subroutine  RKUTTA 1 
Subroutine  DERIV 


Subroutine  DERIV1 

Subroutine  FINISH 

Subroutine  OUTPUT 

Subroutine  BESSEL 
Subroutine  SKETCH 

Subroutine  APPEND 

Subroutine  NOW 


Provides  X  -   I  equation  for  integration  for  a 
particular  potential. 

Computes  phase  shifts,  cross  sections,  and 
normalized  cross  sections. 

Prints  out  resulting  r0/a  ratio,  V0,  rn,  energy, 
phase  shifts,  and  cross  sections. 
Calculates  bessel  functions 
Prepares  for  and  calls  DRAW  to  plot  all  cross 
sections  on  one  graph  for  comparison. 
Stores  base  cross  section  for  normalization 
in  FINISH. 

Prepares  for  and  calls  DRAW  to  plot  all  norma- 
lized cross  sections  on  one  graph  for  comparison, 
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APPENDIX  II 
PROGRAM  DICTIONARY 


EINT 

EINC 

II 

ROS 

AS 

CHANGE 

MANY 
REND 

NUMBER 
MODE 


III 
MOMENTA 


FK 

FK2 

FMANY 

RA 

S 


Initial  value  of  center  of  mass  energy  in  MeV. 

Incremental  value  of  center  of  mass  energy  in  MeV. 

Maximum  number  of  center  of  mass  energy  values. 

Base  effective  range  In  cm. 

Base  scattering  length  In  cm. 

Percentage  Incremental  change  in  base  effective  range 

or  base  scattering  length  in  decimal  form. 

Number  of  +  or  -  incremental  charges  to  be  made. 

Radius  at  which  integrated  solution  is  matched  to 

exact  solution  for  phase  shift  determination  in  cm. 

Number  of  integration  steps. 

Print  out  control  for  program  monitoring 

■  0  for  result  print  out  only 
z   1  for  minor  progress  report 

■  2  for  print  out  of  all  calculations 
Energy  subscript  for  wave  function  print  out. 
z   0   5U  0  calculations  only 

si   |s  1  calculations  only 

*  2  both   |:  0  and   £s  1  calculations. 

Wave  number,  k 

2 

Wave  number  squared,  k  . 

MANY  in  floating  point. 

Ratio  of  effective  range  to  scattering  length,  rQ/a. 

Well  depth  parameter,  s 
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B 

POT 

RN 

ARGB 

VALUE 

VALJO 

VALJ1 

VALJ2 

VALNO 

VALN1 

VALN2 

A 

V 

FNUMBER 

STEP 

R 

YVARS(I) 

WARS  (2) 

XVAR 

Y 

YDOT 

COMPARO 

COMPAR1 

XC 


Intrinsic  range,  b 

Well  depth,  V0 

Well  range,  rn 

Bessel  function  argument,  kr 

Output  array  of  bessel  functions 

Spherical  j|  s  0  bessel  function,  j 

Spherical  j[  s  1  bessel  function,  j| 

Spherical  (s  2  bessel  function,  J£ 

Spherical  j ■  0  neumman  function,  Hq 

Spherical  §  s  1  neumann  function,  n. 

Spherical  J|  ■  2  neumann  function,  n£ 

A  particular  potential  range  from  array,  RN 

A  particular  potential  depth  from  array,  POT 

Number  of  integration  steps  in  floating  point 

Integration  step  size 

Array  of  radius  values  at  each  integration  step 

Wave  function 

Slope  of  wave  functions 

Radius  variable  of  integration 

Wave  function  at  end  value  of  integration  step  for  (:  0. 

Slope  of  wave  function  at  end  value  of  integration 

step  for  II  ■  0. 

Logrithmic  derivative  of  wave  function  at  end  of 

integration  range  for  f  s  0. 

Logrithmic  derivative  of  wave  function  at  end  of 

integration  range  for  0-  1. 

Runge  kutta  weighted  radius  step 
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YC 

DY 

YL 

YLDOT 

DELTAO 

DELTA  1 

XSECO 

XSEC1 

TOTAL 

BETTER 

FNORM 

NUMBERR 


Dependent  variable  for  XC 

Derivative  of  dependent  variable  for  XC 

Same  as  Y  for  J(s  1. 

Same  as  YDOT  for  jj  s  1  . 

J  s  0  phase  shift. 

I    s  1  phase  shift. 

A  s  0  cross  section 

\   »  1  cross  section 
Sum  of  singlet  and  triplet  cross  sections 
Normalized  $=0  cross  sections 
Base  cross  section  to  be  used  for  normalization 
Number  of  simultaneous  first  order  differential  equa< 
tions  to  be  integrated. 
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APPENDIX  III 
PROGRAM  LISTING 


6+ 


MAIN  PROGRAM 

DIMENSION  E(100) »FK(100) »FK2( 100) ,RN( 100) ,VALJ0( 100) .VALJK 100) 
1,VALJ2( 100) ♦  VALNO(IOO),  VALNl ( 100 ) » VALN2 ( 100 ) »  R(IOOO) 


Y( 1000) , 

2YDOTU000),  COMPAR0(  100)  »  C0MPAR1  (  100  )  ,  YL  (  1  000  )  ,YLDOT  (  1000  ),  POT  (  1 
300) ,DELTAO( 100 ) ,DELT Al ( 100 ) ,XSEC0(100) ,XSEC1 (100) ♦ TOTAL ( 100) »RA(10 
40)  ,FNORM( 100), BETTER (10,100) 

COMMON  E,EINT,EINC,  FK»  FK2  ,  RINT,  RINC,  RN,  A,  II,  LL ,  NN,  MODE, 
1REND,VALJ0,VALJ1,VALJ2,VALN0,VALN1,VALN2,JJ,M,  NUMBER,  R,I ,L,V,K, 
2N,Y,  YDOT,  YL,  YLDOT,  III,  COMPARO,  COMPAR 1 ,POT , V  I  NT , VI NC .MOMENTA , 
3DELTAO,DELTA1,XSECO,XSEC1,TOTAL,RA,ROS,AS,CHANGE,MANY,S,B,FNORM,BE 
4TTER 

CALL  INPUT2 

CALL  READY2 

A=RN(MANY+1) 

V=P0T(MANY+l)/.41499E-24 

CALL  FIGURE 

CALL  FINISH 

CALL  APPEND 

DO  1  L=1,LL 

A=RN(L) 

V=  POT(L) /.41499E-24 

CALL  FIGURE  , 

CALL  FINISH 

CALL  0UTPUT2 

CALL  SKETCH(E,XSEC0,MOD,NUMPTS,II ,L,LL) 
.  CONTINUE 

CALL  NPW(E, BETTER, II, L»LL) 

END 


SUBROUTINE  INPUT2 

DIMENSION  E(  100)  ,FK(  100)  ,FK2(100)  ,RN( 100)  ,VALJ0(100)  .VALJK100) 
1 ,VALJ2( 100)  ,  VALNO(IOO),  VALN1 ( 100 ) ♦ VALN2 ( 100 ) »  R(1000)»  Y(1000)t 
2YDOK1000),  COMPARO( 100) ♦  COMPARl ( 100 ) »  YL ( 1000 ) ,YLDOT ( 1000 ) »POT < 1 
300) »DELTA0( 100) ,DELTA1 ( 100 ) ♦XSECO ( 100 ) ,XSEC1<100) •TOTAL < 100) »RA(10 
40) ,FNORM( 100) » BETTER ( 10»100) 

COMMON  E,EINT,EINC,  FK,  FK2.  RINT.  RINC,  RN,  A.  II,  LL,  NN,  MODE* 
1REND»VALJ0,VALJ1,VALJ2»VALN0,VALN1,VALN2»JJ»M,  NUMBER.-  R,I,L,V,K. 
2:N»Y»  YDOT,  YL,  YLDOT  ♦  III*  COMPARO,  COMPARl ,  POT  ,V  INT  »VINC  »MOMENTAt 
3DELTA0,DELTA1,XSEC0,XSEC1» TOTAL, RA ,ROS , AS, CHANGE, MANY ,S.B,FNORM, BE 
4TTER 

READ  1  ,EINT,EINC,ir 

1  FORMAT(2E20. 4,110) 

READ  2, ROS, AS, CHANGE, MANY 

2  ;FORMAT(2E20.6,F10.4,  110)       ^ 
READ  3,REND, NUMBER 

3  FORMAT(E20.4,I10) 
READ  4, MODE 

4  FORMAT( 110) 
READ  5,111 

5  FORMAT ( 110) 
READ  6, MOMENTA 

6  FORMAT(IIO) 
IF(MODE  -  1)  8,8,7 

7  PRINT  1,EINT,EINC,II 
PRINT  2,R0S, AS, CHANGE, MANY 
PRINT  3, REND, NUMBER; 
PRINT  4, MODE 

PRINT  5,111 
PRINT  6, MOMENTA 

8  RETURN 
END 


I 
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SUBROUTINE  FIGURE 

DIMENSION  E(100).FK(  100)  ,FK2d00)  ,RN(100)  ,VALJ0(100)  ♦VALJK100) 
1  ,VALJ2( 100)  ,  VALNO(IOO),  VALN1 ( 100 ) ,VALN2 ( 100 ) »  R(IOOO),  Y(1000)» 
2YDOT(1000),  COMPAR0(  100)»  COMPARl ( 100 ) ,  YL ( 1000 ) ,YLDOT ( 1000 ) ,POT { 1 
300) ,DELTA0( 100) ,DELTA1 ( 100  ), XSECO ( 100 ) ,XSEC1(100) tTOTAL(lOO) ,RA(10 
40) ,FNORM(10  0) » BETTER ( 10 » 100 ) 

DIMENSION  YVARSOO) 

COMMON  E,EINT»EINC»  FK»  FK2»  RINT,  RINC,  RN,  A»  II,  LL,  NN,  MODE, 
1REND»VALJ0,VALJ1,VALJ2»VALN0»VALN1,VALN2»JJ»M,  NUMBER t  R.I »L,VtK, 
2N,Y,  YDOT,  YL»  YLDOT  ♦  III,  COMPARO ,  COMPARl , POT , VI  NT ,VINC .MOMENTA, 
3DELT AO, DELT Al, XSECO, XS EC  1, TOT AL,RA,ROS, AS, CHANGE, MAN Y,S,B,FNORM, BE 
4TTER 

;FNUMBER  =  NUMBER 

STEP  =  REND  /  FNUMBER 

R{ 1)  =  0. 

NN  =  NUMBER  +  1  r 

DO  299  N=2,NN 
299  R(N)  =  R(N-l)  +  STEP 

DO  292  K=2,II 
'  YVARS(l)  =  0, 

YVARS(2)  =  1. 

IF(MOMENTA  -  1)  283,284,283 

283  CALL  RKUTTA(2,XVAR,YVARS,STEP) 
IF(K-III)  293,297,293 

297  PRINT  296 

296  FORMAT  (13X,7H  ENERGY , 13X , 7H  RADIUS, 6X , 14H  WAVE  FUNCTION) 

PRINT  295, E(K) 
295  FORMAT(F20.2) 

PRINT  294,  (R(N),  Y.(N),  N=1,NN) 
294  FORMAT  ( E40.4,E20.4 ) 
293  COMPARO ( K ) =YDOT (NN-1 )/Y(NN-l) 

284  IF  (MOMENTA  -  1)  2^2,291,291 

291  YVARS(l)  =  0. 
YVARS(2)  =  1. 

CALL-RKUTTA1(2,XVAR,YVARS,STEP) 
IF(K-III)  285,289,285 

289  PRINT  288 

288  FORMAT  ( 13X ,7H  ENERGY , 13X,7H  RADIUS, 6X, 18H  WAVE  FUNCTION  L=l)  " 

PRINT  287, E(K) 
287  FORMAT(F20.2) 

PRINT  286,  (R(N) ,YL(N)  ,N=1,NN) 
286  FORMAT(E40.4,E20.4) 

28  5  COMPARl  (  K)=YLDOT(  NN-1  )/YL(  NN-1). 

292  CONTINUE 
RETURN 
END 
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SUBROUTINE       RKUTTA ( NUMBERR ,XVAR , WARS. STEP ) 

DIMENSION    E( 100) ,FK( 100) ,FK2Q00) »RN( 100) .VALJ0Q00) .VALJK 100) 
1,VALJ2( 100) ,    VALNO(IOO),    VALNl ( 100 ) . VALN2 ( 100 ) ♦    R(IOOO),    Y(IOOO). 
2YDOTU00O),    COMPAROC 100) ,    COMPARl ( 100 ) »    YL ( 1000 ) .YLDOT ( 1000 )  .POT ( 1 
300)  .DELTAOC  100 ) .DELTAl < 100 )  ,XSEC0<100)  .XSEC1U00)  .TOTAL  (100)  tRAdO 
40)  ,FNORM( 100) , BETTER ( 10.100) 

DIMENSION    YVARSOO)  »  AK  (  4,  30  )  .  DY  (  30  )  .  YC  (  30  )  .C  (  4  ) 

COMMON    E.EINT.EINC.    FK,    FK2.    RINT.    RINC.    RN ,    A.     II,    LL ♦    NN ,    MODE. 
1REND.VALJ0,VALJ1,VALJ2.VALN0.VALN1.VALN2»JJ»M,    NUMBER.    R.I .L.V.K, 
2N,Y,    YDOT,    YL,    YLDOT  ,     III.    COMPARO,    COMPARl , POT ,Vl NT , VI NC, MOMENTA , 
3DELTAO,DELTA1,XSECO,XSEC1»TOTAL,RA,ROS.AS»CHANGE.MANY.S,B.FNORM.BE 
4TTER 

C( 1 )=0.0 

C(2)=0.5 

C(3)=0.5 

C(4)=1.0 

DO    298    N=1,NN 

XVAR       =    R(N) 

DO    1    1=1.4 

XC=XVAR    +    C( I )*STEP 

DO    2    J=l, NUMBERR 

2  YC(J)    =YVARS(J)    +    C( I )*AK( I-l.J) 
CALL    DERIV(XCYC.DY) 

DO    1    J=1.NUM8ERR 
1    AK( I »J)    =    STEP    *    DY( J) 
DO    3    J    =    1, NUMBERR 

3  YVARS(J)     =YVARS(J)     +     ( AK ( 1 , J ) +2.*AK ( 2 » J ) +2.*AK ( 3 » J ) +AK( 4. J ) ) /6. 
YDOT(N)     =    YVARS<2) 

298    Y(N)     =    YVARS(l) 
RETURN 
END 
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SUBROUTINE  RKUTTA1 ( NUMBERR , XVAR ,YVARS » STEP ) 

DIMENSION  E( 100) ,FK( 100) ,FK2(100) ,RN(100) ,VALJ0<100) ,VALJ1( 100) 
1»VALJ2(  100) ,  VALNO(IOO),  VALN1 ( 100 )  ,VALN2 ( 100 ) »  R(IOOO),  Y(1000)» 
2YDOTU000),  COMPAROf 100) »  COMPARl ( 100) .  YL ( 1000 ) ,YLDOT ( 1000 ) »POT (1 
3  00) ,DELTAO( 100 )  ,DELT Al ( 100  ) ,XSEC0(100) tXSECK 100) tTOTAL(lOO) »RA(1C 
40) »FNORM( 10  0) » BETTER ( 10 , 100 ) 

'DIMENSION  YVARSOO)  ,  AK  (  4.30  )  »  DY  (  30  )  ,  YC  (  30  )  ,C  (  4  ) 

COMMON  E,EINT,EINC,  FK,  FK2»  RINT.  RlNCt  RN »  At  II,  LL»  NN ♦  MODE, 
1REND,VALJ0,VALJ1,VALJ2,VALN0,VALN1,VALN2,JJ,M,  NUMBER,  R,I ,L,V,K, 
2N,Y,  YDOT,  YL,  YLDOT  ,  III,  COMPARO ,  COMPARl , POT ,V INT ,VINC , MOMENTA, 
3DELTAO,DELTA1,XSECO,XSEC1,TOTAL,RA,ROS,AS,CHANGE,MANY,S,B,FNORM,B£ 
4TTER 

C( 1)=0.0 

C(2)=0.5 

C(3)=0.5 

C(4)=1.0 

DO  290  N=1,NN 

XVAR  =  R(N) 
I  DO  1  1  =  1,4 
;XC=XVAR  +  C( I )*STEP 
;DO  2  J=l, NUMBERR   . 

2  YC(J)  =YVARS(J)  +  C< I )*AK< I-l.J) 
CALL  DERIVKXCYCDY) 

DO  1  J=l, NUMBERR 
1  AK( I, J)  =  STEP  *  DY( J) 
DO  3  J  =  1, NUMBERR' 

3  YVARS(J)  =YVARS(J)  +  ( AK ( 1 , J ) +2.*AK ( 2 , J )+2.*AK ( 3 , J ) +AK( 4, J ) ) /6 • 
YL(N)  =  YVARS(l) 

290jYLDOT(N)  =  YVARS(2) 
RETURN 
END 
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SUBROUTINE  FINISH 

DIMENSION  E(  100)  ,FK(  100)  ,FK2(100)  ,RN( 100)  ,VALJ0(100)  tVALJH  100) 
1  ,VALJ2( 100)  ,  VALNO(IOO),  VALN1 ( 100 )  , VALN2 ( 100 ) »  R(1000)»  Y(1000)» 
2YDOTU000),  COMPAR0(  100)  »  COMPARl  ( 100  )  ♦  YL ( 1000 ) ,YLD0T ( 1000 ) »P0T ( 1 
3  00) »DELTA0( 100 ) ,DELT Al ( ICO ) , XSECO ( 100) ,XS£C1 (100) .TOTAL < 100) iRAIlO 
40) , FN0RM( 100) , BETTER ( 10*100) 

COMMON  E»EINT,EINC,  FK,  FK2,  RINT,  RINC,  RN,  A*  II,  LL,  NN ,  MODE, 
1REND,VALJ0,VALJ1,VALJ2,VALN0»VALN1»VALN2»JJ,M,  NUMBER,  R,I ,L,V,K, 
2N,Y,  YDOT,  YL,  YLDOT ,  III,  COMPARO,  COMPARl , POT ,V I  NT , VI NC .MOMENTA , 
3 DELT A0,DELTA1,XSEC0,XSEC1, TOTAL, RA,ROS, AS, CHANGE, MAN Y,S,B,FNORM, BE 
4TTER 
IF(MOMENTA  -  1)  399,397,399 
399  DO  398  1=2,11 

DELTAO( I )=ATANF( ( REND* (COMPARO ( I)*VALJO( I)+FK( I )*VALJ1 ( I ) )- 
1VALJ0( I ) ) /( REND* ( COMPARO ( I )*VALNO( I )+FK( I )*VAlN1( I ) )- 
2VALN0( I ) ) ) 
398  XSECO( I )=4.*3.1416/FK2( I ) *SI NF ( DELTAO ( I ) )**2 
397  IF(M0MENTA-1)394,396,396 
396  DO  395  1=2,11 

DELTAK  I)=ATANF( (REND*(FK( I )/3.*VALJ0( I )~2./3.*FK(  I  )*  • 
1VALJ2( I )-COMPARl( I )*VALJ1( I ) )+VALJl( I ) ) /(REND*(FK( I )/3.* 
2VALN0( I )-2./3.*FK( I ) *VALN2 ( I ) -COMPARl U )*VALN1(I ) )+VALNl( I ) ) ) 
39  5  XSECK I )=12.*3.1416/FK2( I ) *SINF ( DELTA1 ( I ) )**2 
390  DO  389  I  =  2,11 

389  TOTAL(I)  =  .25  *  XSECO  (  I  )  +  ,75  *  XSECKI) 
394  DO  900  1=2,11       J 

BETTER(L,I)  =  XSECO ( I ) /FNORM ( I ) 
900  CONTINUE  j. 

RETURN 
END 


67 


SUBROUTINE  OUTPUT2 
! pi  MENS  I  ON  E(100) ,FK( 10  0) »FK2<100) t RN ( 100 )»VALJO < 100 )  »VALJ1(100) 
l'»VALJ2(100)  »  VALNO(IOO),  VALN1  (  100  )  ,  VALN2  (  100  )  »  R(1000)»  Y(IOOO), 
2YDOT(1000),  COMPAR0(  100).,  COMPARl ( 100 ) ♦  YL ( 1000 ) »YLDOT ( 1000 ) tPOT ( 1 
300)  ,DELTA0(  100 )  ,DELTA1  ( 100  )  tXSECOUOO)  ,XSEC1  (100)  .TOTAL  (100)  ,RA(10 
40) ♦FNORM( 100) » BETTER ( 10, 100 ) 

COMMON  E,EINT,EINC,  FK,  FK2,  RINT,  RINC,  RN,  A,  II,  LL ,  NN,  MODE, 
1REND,VALJ0,VALJ1,VALJ2,VALN0,VALN1,VALN2,JJ,M»  NUMBER,  R»I»L»V,K, 
2N,Y,  YDOT,  YL,  YLDOT,  III,  COMPARO,  COMPARl , POT , VI  NT ,VINC .MOMENTA, 
3DELT AO, DELT Al » XSECO » XSEC1, TOT AL,RA,ROS, AS, CHANGE, MAN Y»S,B,FNORM, BE 
4TTER 
PRINT  600,RA(L) 
600  F0RMATO4H  EFFECTIVE  RANGE/SCATTERING  LENGTH  ,/,  E34.  10  ) 

PRINT  599,A,POT(L) 
599  FORMAT  (13X,7H  RADIUS ,4X , 16H  POTENTIAL  DEPTH ,/ ,E20.4 ,E20«4 ) 

PRINT  598 
598  F0RMAT(7H  ENERGY, 12H  PHASE  ANGLE, 15H  CROSS  SECT  ION, 12H'  PHASE  ANGLE 
1,15H  CROSS  SECTION, 2X,6H  TOTAL ,/, 1 IX ,4H  L=0,7X ,4H  L=0» 10X, 4H  L=l, 
28X,4H  L=1,8X,6H  CROSS,/ ,63X,8H  SECTION) 

PRINT  597, (E( I ) ,DELTAO( I ) ,XSECO( I ) ,DELTA1( I) ,XSEC1(I ) .TOTAL ( I) ,1=2 
1,11  ) 
597  F0RMAT(F7.2»E12.4,E15.4,E12.4,E15.4,E11.4) 

PRINT  596 
596  F0RMAT(7H  ENERGY , 5X, 25H  NORMALIZED  CROSS  SECTION) 

PRINT  595, (E( I ) ,BETTER(L,I ) ,1=2,11 )   -• 
595  FORMAT(F7.2,E20.8) 
RETURN 
END 
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SUBROUTINE  BESSEL  ( X »N »VALUE ♦ I NDEX ,MM ) 
DIMENSION  VALUE (40) , BJ(lOO)  »F( 100) 
IF  (X)  10,10,9 
9  IK=X/2.0 
K=N+20+IK 
IF(K-IOO)  300,299,299 

299  PRINT  298 

298  F0RMAT(66H  ARGUMENT  PLUS  ORDER  TOO  LARGE.  INCREASE  DIMENSIONS  AND 
1RUN  AGAIN.) 
RETURN 

300  CONTINUE 
Z=1.0/X 

PI=3. 141592653 

GOTO (40, 40, 60, 60, 80, 90, 100, 110) , I NDEX 

10  PRINT  8,X 

80FORMAT  ( IX, 63HSUBROUT INE  BESSEL  DOES  NOT  WORK  FOR  X=0  OR  LESS  THAN 
1;  0,  HERE  X=  E15.8) 
RETURN 

40  BJ(K+1)=0.0 
BJ(K)=10.0E-30 
L  =  K 

11  FL=L 
BJ(L-1)=2.0*FL*Z*BJ(L)-BJ(L+1) 
IF(L-l)  12,12,13 

;  13  L  =  L-1 

GO  TO  11 

12  SUM=0.0 
IK  =  K 

DO  14  I=2,IK,2 

14  SUM=SUM+BJ( I  ) 
C=1.0/(BJ(0)+2.0*SUM) 
GO  TO(42,43) , INDEX  , 

43  DO  45  1=0, K 
45  BJ( I )=C*BJ( I  ) 

GO  TO  54 
42  DO  41  1=0, N 

41  VALUE( I )=C*BJ( I ) 
RETURN 

50  GO  TO  40 

54  EC=. 5772156649 

SUM=0.0 

IK=K/2 

DO    15    1=1,  IK 

FI  =  I 

15  SUM=SUM+( (-1.0)**( I-1)*BJ(2*I ) )/FI 
F(0)=(2.0/PI)*(BJ(0)*(LOGF(X/2.0)+EC)+2.0*SUM) 
F(l )=(BJ(l)*F(0)-2.0/(PI*X) )/BJ(0) 

IN=N-1 

DO    16    1=1, IN 

,FI  =  I 
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16  F( I+1)=2.0*FI*Z*F( I )-F( 1-1 ) 
DO  51  1=0, N 

51  VALUE( I )=F( I ) 
RETURN 

60  ;  BJ(K+1)=0. 

BJ(K)*10.E-30 

L  =  K 
19  FL=L 

BJ ( L-l ) =2 .*FL*Z*BJ ( L ) +BJ ( L+l ) 

IF(L-l)17f 17,18 
18  L=L-1 

GO  TO  19 

17  SUM=0.0 

DO  21  1=1, K 

21  SUM=SUM+BJ( I ) 
C=(EXPF(X) )/(BJ(0)+2.*SUM) 
GO  TO( 1,1,61,62) , INDEX 

1  PRINT  1000 
RETURN 
1000  FORMAT(40H  CARD  MISSING.  IGNORE  RESULTS 
62,  DO  64  1=0,1 

64  BJ( I )=C*BJ( I ) 
i  ;GO  TO  71 

61  DO  6  5  I=0,N 

65  VALUE( I )=C*BJ( I ) 
RETURN 

.  70  GO  TO  60 
71  IF(5.0-X)  120,121,121 

121  M=10 

GO  TO  122 
120  M=5 

122  F  (0)=GAUSS(0.0,1.0,M,X) 
F(1)=(1./X-F(0)*BJ(1) )/BJ(0) 
IN=N-1  , 

DO  22  1=1, IN 
FI  =  I 

22  F  ( I+1)=2.0*FI*Z*F  (I)+F  (1-1) 
DO  72  I=0,N 

72,  VALUE(  I  )  =  F  (  I  ) 

RETURN 
80  F(K+1) =0.0 

F  (K)=10.0E-30 

L  =  K 
25  FL=L 

F  (L-1)=(2.0*FL+1.0)*Z*F  ( L) -F  (L+l) 

IF  (L-l)   23,23,24 
24  L=L-1 

GO  TO  25 

23  C=SINF(X)/(X*F  (0) ) 
IF(MM)200,200,201 

201  SQ=SQRTF(2.*X/PI ) 
DO  202  1=0, N 

202  VALUE( I )=C*SO*F  (I) 
RETURN 
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200  DO  26  1=0, N 

26  VALUEf  I  )=C*F  (  I  ) 
RETURN 

90  F  (0)=-COSF(X)/X 

F  ( 1)=-SINF(X)/X-C0SF(X)/X**2 

IN=N-1 

DO  27  1=1, IN 

FI  =  I 

27  F(I+1)=(2.*FI+1.0)*Z*F(I)-F(I-1) 
IF(MM)  203,203,204 

204  SQ=SQRTF(2.*X/PI ) 
DO  205  I=0,N 

205  VALUE(I)=F  ( I )*SQ*<-1 . )**<-I-l ) 
RETURN 

203  DO  91  1=0, N 

91  VALUEf  I )=F  (  I  ) 
RETURN 

100  F  (K+1)=0.0 

F  (K)=10.0E-30 
L  =  K 

30  FL=L 

F  (L-1)=(2.0*FL+1.0)*Z*F  (L)+F  (L+l) 
IF  (L-l)  28,28,29 
29  L=L-1 

GO  TO  30 

28  C=SINH  (X)/(X*F  (0) ) 
IF(MM)  206,206,207 

207  SQ=SQRTF(2.*X/PI) 
DO  208  I=0,N 

208  VALUEf I )=C*SQ*F  (  I  ) 
RETURN 

206  DO  101  1=0, N 

101  VALUEf I )=C*F  ( I ) 
RETURN 

110  JF  (0)=-PI*EXPF(-X)/(2.0*X) 

jF  (1  )  =  (PI*(X  +  1.0)*EXPF(-X)  )/(2.0*X**2) 

TN=N-1 

DO  31  1=1, IN 

FI  =  I 

31  F  (I+1)=F  ( I-1)-(2.0*FI+1.0)*Z*F  (I) 
IF(MM)  209,209,210 

210  SQ  =  SQRTF(2.*X/PI  ) 
DO  211  1=0, N 

211  VALUE(.I)=F  (  I  )*SQ*(-1.  )**(-I-l) 
RETURN 

209  DO  111  1=0, N 

111  VALUE( I )=F  { I ) 
RETURN 

!END 

FUNCTION  GAUSS   NUMERICAL  INTERGRATION 

FUNCTION  GAUSS  REQUIRES  FUNCTION  GAUSSF(X)  NUMERICAL  INTERGRATION 
FUNCTION  GAUSS(G,B,M,X) 
.DIMENSION  R(10) ,U(10) 
,U(1)=+. 0744371695 
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123 


;U( 
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U( 
U( 
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D  = 
C  = 
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RE 
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)/FN 
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)*GAUSSF< (Y- 


■A)*U( I )+(A+Y)/2.0tX) )*(Y-A) 
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1  J' 
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2  I 
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D+C 
Y 

USS=GAUSS  +C 
TURN 
D 
NCTT 

(1.0, 

USSF' 

TURN 

USSF=  (EXPF(-X*COSH(1.0/U-1.0) ) ) /U**2 

D 


ON  GAUSSF(U»X) 
i/U-710. 0)123, 124, 124 
=  0.0 
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SUBROUTINE  SKETCH ( E ♦ XSECO .MOD.NUMPTS » I  I »L,LL) 
DIMENSION  E( 100) ,XSEC0(100)  ,XSEC1(100) . TOTAL (100) 
DIMENSION  X(900),Y(9  00)»ITITLE(12) 
LABEL  =  4H 
DO  1  1*1*12 

1  ITITLE( I)  =  8H 
ITITLE( 1)=8H0SG00D»F 
ITITLE(2)=8H.B, ,B0X 
ITITLE(3)=8H09 
ITITLE(A)=8HEXP0NENT 
ITITLE(5)=8HIAL  WELL 
ITITLE(6)=8H  RUN  15 
ITITLE(7)=8H  GRAPH  0 
ITITLE(8)=8HF  CROSS 
ITITLE(9)=8HSECTI0N 
ITITLE( 10)=8HVS  ENERG 
ITITLE( 11)=8HY  FOR  L= 
ITITLE(12)  =  8H0 
NUMPTS  =11-1 
IF(L-1)2»2»3 

2  MOD  =  1 
GO  TO  6 

3  IF(LL-L)  5»5»4 

4  MOD  =  2 
GO  TO  6 

5  MOD  =  3 

6  CONTINUE 

CALL  DRAW <NUMPTS>E .XSECO » MOD* 0> LABEL » I TI TLE » • 3 »0 ,0 »0 »0 »0»9 . 15 » 1 »L 
1AST) 
RETURN 
END 
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SUBROUTINE  APPEND 

DIMENSION  E( 100) ,FK( 100) ,FK2(100) ,RN(100) .VALJOUOO) »VALJ1( 100) 
1 »VALJ2( 100) ,  VALNO(IOO),  VALNl ( 100 ) ,VALN2 { 100 ) »  R(1000)»  Y(1000)» 
2YDOT(1000),  COMPAR0( 100) ,  COMPARl ( 100 ) ,  YL ( 1000 ) ,YLDOT ( 1000 )  ,POT ( 1 
3  00) ,DELTAO( 100 ) ,DELT Al ( 100 ) ,XSEC0(100) , XSECl (100) » TOTAL (100) ,RA(10 
40) ,FNORM( 100) » BETTER (10,100) 

COMMON  E,EINT,EINC,  FK,  FK2»  RINT,  RINC»  RN,  A»  II,  LL,  NN»  MODE, 
lREND,VALJ0,VALJl,VALJ2,VALN0,VALNl,VALN2,JJ,M,  NUMBER,  R,I,L,V,K, 
2N,Y,  YDOT,  YL,  YLDOT,  III,  COMPARO,  COMPARl , POT , VI  NT , VI NC , MOMENTA, 
3 DELT AO, DELTA1,XSEC0, XSECl, TOT AL,RA,ROS, AS, CHANGE, MAN Y,S,B,FNORM, BE 
4TTER 

DO  1  1=2,11 
1  FNORM(  I )  =  XSECO(  I  ) 

RETURN 

END 
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SUBROUTINE  NOW( E .BETTER ♦ 1 1 »L , LL) 

DIMENSION  E(100),  BETTER < 1 0 , 100 ) , X ( 900 } , Y ( 900 ) , I T I TLE ( 12 ) 

LABEL  =  4H 

DO  1  1=1,12 
1  ITITLE(I)=8H 

ITITLE( 1)=8H0SG00D,F 

ITITLE(2)=8H.B. .BOX 

ITITLEt 3)=8H09 

ITITLE(7)=8HNORMALIZ 

ITITLE(8)=8HED  CROSS 

ITITLE(9)=8H  SECTION 

ITITLE(10)=8H  FOR  L=0 

NUMPTS  =  I  I  —  1 

DO  7  L=1,LL 

IF(L-l)  2,2,3 
2, M0D=1 

GO  TO  6 

3  IF(LL-L)  5,5,4 

4  MOD  =  2 
GO  TO  6 

5  MOD  =  3 

6  CONTINUE 

DO  20  1=2,11 
20  BETU  )=BETTER(L,I  ) 

CALL  DR AW ( NUMPTS, E, BET,    MOD,0 , LABEL,  IT  I TLE , . 3 ,0 ,0 » 0 ,0 ,0 ,9  ,  15  , 1 ,  L 
1AST) 

7  CONTINUE 
RETURN 

END  :, 
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THE  FOLLOWING  SUBROUTINES  ARE  DEPENDENT  UPON  A  PARTICULAR  POTENTIAL  SHAPE. 
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GAUSSIAN  POTENTIAL 

SUBROUTINE  READY2 

DIMENSION  E(100)»FK(  100)  »FK2(100)  ,RN(100)  ,VALJ0(100)  >VAL_J1  (100) 
1  ,VALJ2( 100)  ,  VALNO(IOO),  VALN1 ( 100 ) , VALN2 ( 100 ) ,  R(1000),  Y(1000), 
2YDOT(1000),  COMPAR0( 100) »  COMPARl ( 100 ) »  YL ( 1 000 ), YLDOT ( 1000 ), POT ( 1 
3  00)  ,DELTAO( 100 )  ,DELT Al ( 100 ) ,XSECO( 100 ) , XSEC1 ( 100 ) , TOTAL { 100)  »RA(10 
40)  ,FNORM( 100)  , BETTER ( 10,100) 

DIMENSION  VALUE(40) 

COMMON  E,EINT,EINC,  FK,  FK2,  RINT,  RINC,  RN,  A»  II,  LL,  NN ,  MODE, 
1REND,VALJ0,VALJ1,VALJ2,VALN0,VALN1,VALN2,JJ*M,  NUMBER,  R,I ,L,V,K, 
2N,Y,  YDOT,  YL,  YLDOT,  III,  COMPARO,  COMPARl , POT ,V I  NT ,V I NC , MOMENTA , 
3DELTAO,DELTA1,XSECO,XSEC1,TOTAL,RA,ROS,AS,CHANGE,MANY»S,B,FNORM,BE 
4TTER 

E( 1 )  =  EINT 

DO  199  I  =  2,11 

E(  I  )  =  E( 1-1 )  +  EINC 

FK(I)  =  SQRTF(E( I ) / .41499E-24 ) 
199  FK2( I )=  FK(  I )  *  FK( I  ) 

LL=2*MANY  +  1 

JJ  =  LL 

FMANY  =  MANY 

GAUSSIAN  DISTRIBUTION 

RA(1)  =  R0S/AS*(1.-  FMANY*CHANGE) 

DO  170  L=1,LL 

S=(1.-2.0  3*RA(L)+.58  7*RA(L)**2)/(1.-2.81*RA(L)+1.805*RA(L)**2) 

ss=s-i. 

B=RA(L)*AS/( l.-.641*SS+.568*SS**2-.721*SS**3+1.172*SS**4-1.015*SS* 
1*5) 
POT(L)=S*.2  2920  8E-23/B**2 
RN(L)=SQRTF(B**2/2.0604) 

170  RA(L+1)  =  RA(L)  +  CHANGE  *ROS/AS 
IF(MODE  -  1)  190,192,197 

197  PRINT  196 

196  FORMAT  (13X,7H  ENERGY , 5X , 1 2H  WAVE  NUMBER , 3X, 20H  WAVE  NUMBER  SQUARE 
ID) 

PRINT  195  ,(E(I),  FK(  I  )  ,FK2(I  )  ,1=2,11) 
195  FORMAT  (3E20.4) 

PRINT  171 

171  F0RMAT(34H  EFFECTIVE  RANGE/SCATTERING  LENGTH) 
PRINT  172,(RA(L)  L=1,LL) 

172  FORMAT(E34.10) 
PRINT  194 

194  FORMAT  (13X,7H  RADIUS) 

PRINT  193,  (RN(L),  L=1.LL) 
193  FORMAT  (E20.4) 
PRINT  1192 
1192  FORMAT  ( 16H  POTENTIAL  DEPTH) 

PRINT  2192, (POT(J) ,J=1,JJ) 
2192  F0RMAT(F16.2) 
192  PRINT  191 
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191  FORMAT  (  34-H  RADIUS  AND  ENERGY  VALUES  COMPUTED) 
190  DO  187  1=2,11 

AR6B      =  FK( I )  *  REND 

CALL  BESSEL(  ARGB , 2 , VALUE » 5 ,  0) 

VALJ0(  I  )  =  VALUE(O) 

VALJK  I  )  -  VALUE*  1  ) 

VALJ2(  I  )  =  VALUE(2) 

CALL  BESSEL  ( ARGB , 2 , VALUE,  6,0) 

VALNO( I )     =    VALUE(O) 

VALNK  I  )     =    VALUE(  1  ) 

VALN2(  I )     =    VALUE(2  ) 

IFfMODE    -    1)     187,187,189 
189    PRINT    188,1 
188    FORMATUM5) 
187    CONTINUE 

IF(MODE    -    1)     184,186,186 
186    PRINT    185 

185    FORMAT (     26H    BESSEL    FUNCTIONS    COMPUTED) 
184    RETURN 

END 
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EXPONENTIAL  POTENTIAL 


SUBROUTINE  READY2 

DIMENSION  E( 100) ,FK( 100) »FK2(100) ♦ RN < 100 ) ♦ VAL JO ( 100 ) »VALJ1 (100) 
1  ,VALJ2( 100)  »  VALN0(100)»  VALN1 ( 100 )  » VALN2 ( 100 )  .  R(1000)t  Y(1000)> 
2YDOT(1000),  COMPAROt  100)  *  COMPAR 1 ( 100 )  ,  YL ( 1 000 ), YLDOT ( 1000 ) *POT ( 1 
300)  ,DELTA0( 100)  »DELTA1 (100) ,XSEC0(100) ,XSEC1( 100) *  TOTAL (100)  »RA(10 
40)  ,FNORM(100) ,BETTER(10>100) 

DIMENSION  VALUE(40) 

COMMON  E,EINT,EINC,  FK  ♦  FK2,  RINT»  RINC,  RN ,  A,  II,  LL »  NN,  MODE, 
1REND,VALJ0,VALJ1»VALJ2»VALN0»VALN1 »VALN2»JJ»M,  NUMBER  ♦  R»I »L,V,K, 
2N,Y,  YDOT,  YL,  YLDOT,  III.  COMPARO,  COMPAR1 , POT , V  I  NT ,VI NC ♦ MOMENTA , 
3DELTAO,DELTA1,XSECO,XSEC1,TOTAL,RA,ROS»AS»CHANGE,MANY»S,B»FNORM»BE 
4TTER 

E(l)  =  EINT 

DO  199  I  =  2,11 

E(I)  =  E(I-l)  +  EINC 

FK(I)  =  SQRTF(E( I > /.41499E-24) 
199  FK2( I )=  FK(  I )  *  FK (  I  ) 

LL=2*MANY  +  1 

JJ  =  LL 

EXPONENTIAL  WELL 

FMANY  =  MANY 

RA(1 )=ROS/AS*( l.+FMANY*CHANGE) 

DO  170  L=1,LL 

S=(1.-2.9739*RA(L)+1.3602*RA(L)**2)/(1.-3.7134*RA(L)+3.1177*RA(L)* 
1*2) 

SS=S-1. 

B=RA(L)*AS/( l.-.904*SS+.745#SS**2-.542*SS**3+.174*SS**4) 

POT(L) =S*.751541E-23/B**2 

RN(L)=B/3.5412 

170  RA(L+1)  =  RA(L)  -  CHANGE  *ROS/AS 
IF(MODE  -  1)  190, 192*197 

197  PRINT  196 

196  FORMAT  (13X.7H  ENERGY , 5X , 12H  WAVE  NUMBER , 3X , 20H  WAVE  NUMBER  SQUARE 
ID) 
■  PRINT  195  ,(E(I),  FK( I )  ,FK2(  I  )  .1=2,  II  ) 
195  FORMAT  (3E20.4) 
PRINT  171 

171  FORMAT(34H  EFFECTIVE  RANGE/SCATTERING  LENGTH) 
PRINT  172.(RA(L)  L=1,LL) 

172  FORMAT(E34.10) 
PRINT  194. 

194  FORMAT  (13X,7H  RADIUS) 

PRINT  193,  (RN(L)»  L=1»LL) 
193  FORMAT  (E20.4) 

PRINT  1192 
192  FORMAT  ( 16H  POTENTIAL  DEPTH) 

PRINT  2192, (POT(J) »J=1»JJ) 
192  FORMAT(F16.2) 
192  PRINT  191 
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191  FORMAT (  34H  RADIUS  AND  ENERGY  VALUES  COMPUTED) 
190  DO  187  1=2,11 

ARGB      =  FK( I )  *  REND 

CALL  BESSEL(  ARGB , 2 .VALUE > 5 >  0) 

VALJ0( I )  =  VALUE(O) 

VALJK  I  )  =  VALUE(  1  ) 

VALJ2( I )  =  VALUE(2) 

CALL  BESSEL  ( ARGB , 2 > VALUE »  6.0) 

VALN0( I )  =  VALUF(O) 

VALNK  I  )     =    VALUEU  ) 

VALN2( I )    =    VALUE(2 ) 

IF(MODE    -    1)     187,187,189 
189    PRINT    188,1 
188    F0RMAT(14I5) 
187    CONTINUE 

IF(MODE    -    1)     184,186,186 
186    PRINT    185 

185    FORMATT   26H    BESSEL    FUNCTIONS    COMPUTED) 
184    RETURN 

END 
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YUKAWA  POTENTIAL 

SUBROUTINE  READY2 

DIMENSION  E( 100) ,FK( 10  0) »FK2( 100) ,RN(100) ,VALJ0(100) ,VALJ1( 100) 
1»VALJ2(100) ,  VALNO(IOO),  VALN1 ( 100 ) ♦ VALN2 ( 100 ) »  R(1000),  Y(1000)t 
2YDOT(1000),  C0MPAR0( 100)  ,  COMPARl ( 100 ) ,  YL ( 1000 ) ,YLDOT ( 1000 ) tPOT ( 1 
3  00)  ,DELTAO(  100)  ,DELTA1  (100)  ,XSEC0(  100)  •XSECK100)  t  TOTAL  (100)  »RA(10 
40) »FNORM( 100) » BETTER ( 1 0 , 100 ) 

DIMENSION  VALUE(40) 

COMMON  E,EINT»EINC,  FK »  FK2,  RINT,  RINC,  RN ,  A.  II,  LL,  NN,  MODE, 
1REND,VALJ0,VALJ1,VALJ2,VALN0»VALN1,VALN2»JJ.M,  NUMBER*  R»I ,L,V,K, 
'2N,:Y,  YDOT,  YL,  YLDOT  ,  III,  COMPARO,  COMPARl ,  POT  ,V  INT  ,V  INC  .MOMENTA  , 
3DELTA0,DELTA1,XSEC0,XSEC1, TOTAL, RA,ROS, AS, CHANGE, MAN Y,S,B,FNORM,BE 
4TTER 

E(l)  =  EINT 

DO  199  I  =  2,11 

E(I)  =  E(I-l)  +  EINC 

FK(I)  =  SQRTF(E(  I  ) /.41499E-24) 
199  FK2( I )=  FK( I )  *  FK ( I ) 

LL=2*MANY  +  1 

JJ  =  LL 

FMANY  =  MANY 

RA(l)=ROS/AS*( l.+FMANY*CHANGE) 
:      YUKAWA  WELL 

DO  170  L=1,LL 

S=(1.-3.5093*RA(L)+2.0  3  54*RA(L)**2 ) / ( 1 .-4. 1348*RA ( L ) +3. 7072 *RA ( L) * 
1*2) 

SS=S-1. 

B=RA(L)*AS/( l.-1.369*SS+1.093*SS**2-1.127*SS**3+1.00  5*SS**4-.324*S 

1S**5) 
POT(L)=S*.147585E-23*2.1196/B**2 
RN(L)=B/2. 1.196 

170  RA(L+1)  =  RA(L)  -  CHANGE  *ROS/AS 
IF (MODE  -  1)  190,192,197 

197  PRINT  196 

196  FORMAT  (13X,7H  ENERGY , 5X , 12H  WAVE  NUMBER , 3X, 20H  WAVE  NUMBER  SQUARE 
ID) 

PRINT  195  ,(E(I),  FK( I  )  ,FK2( I )  ,1=2,  II  ) 
195  FORMAT  (3E20.4) 

PRINT  171 

171  FORMAT(34H  EFFECTIVE  RANGE/SCATTERING  LENGTH) 
PRINT  172,(RA(L)  L=1,LL) 

172  FORMAT(E34.10) 
PRINT  194 

194  FORMAT  (13X,7H  RADIUS) 

PRINT  193,  (RN(L) ,  L=1,LL) 
. 193  FORMAT  (E20.4) 

PRINT  1192 
1192  FORMAT  ( 16H  POTENTIAL  DEPTH) 

PRINT  2192* (POT(J) ,J=1,JJ) 
2192  FORMAT(F16.2) 
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192  PRINT  191 

191  FORMAT (  34H  RADIUS  AND  ENERGY  VALUES  COMPUTED) 

190  DO  187  1=2,11 

ARGB      =  FK( I )  *  REND 

CALL  BESSEL(  ARGB , 2 , VALUE, 5 ♦  0) 

VALJO< I )  =  VALUE(O) 

VALJK  I  )  =  VALUE(  1  ) 

VALJ2( I )  =  VALUE<2) 

CALL  BESSEL  ( ARGB» 2 , VALUE.  6,0) 

VALNO( I )  =  VALUE(O) 

VALNK  I  )  =  VALUE(  1  ) 

VALN2( I )  =  VALUE(2 ) 

IF(MODE    -    1)     187,187,189 
189    PRINT    183,1 
188    F0RMATU4I5) 
187    CONTINUE 

IF(MODE  -  1)  184,186,186 
186  PRINT  185 

185  FORMAT*  26H  BESSEL  FUNCTIONS  COMPUTED) 
184  RETURN 

END 
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SQUARE  WELL  POTENTIAL 

SUBROUTINE  READY2 

DIMENSION  E(  100)  ,FK(  100)  ,FK2(l00)  ,  RN  (  100  )  » VAL  JO  (  100  )  »VALJ1(  100) 

1,VALJ2< 100) ,  VALNO(IOO),  VALNl ( 100 )  , VALN2 ( 100  )  ,  R(1000),  Y(1000). 

2YDOT<1000),  COMPAR0( 100) ,  C0MPAR1 ( 100 ) ,  YL ( 1 000 ), YLDOT ( 1000 ) »POT ( 1 

3  00)  ,DELTAO(  100 )  ,DELT  Al  (  100  )  ,XSECO(  100)  ,XSEC1  (100)  ,  TOTAL  (100)  ,RA<10 

40) »FNORM( 100) »BETTER (10,100) 

DIMENSION  VALUE(40) 

COMMON  E,EINT,EINC,  FK»  FK2,  RINT,  RINC,  RN,  A,  II,  LL,  NN,  MODE, 
1REND,VALJ0,VALJ1,VALJ2,VALN0,VALN1,VALN2,JJ,M,  NUMBER,  R, I ,L,V,K, 
2N,Y,  YDOT,  YL,  YLDOT,  III,  COMPARO,  C0MPAR1 , POT , V  I  NT , V  I NC , MOMENTA , 
3DELTAO,DELTA1,XSECO,XSEC1,TOTAL,RA,ROS,AS,CHANGE,MANY,S,B,FNORM,BE 

4TTER 
E( 1)  =  EINT 
DO  199  I  =  2,11 
E( I )  =  E( 1-1  )  +  EINC 
FKU)  =  SQRTF(E( I ) /.41499E-24) 
199  FK2( I )=  FK( I )  *  FK ( I ) 
LL=2*MANY  +  1 
JJ  =  LL 

SQUARE  WELL 
FMANY  =  MANY 

CHANGE  SIGN  TO  -  FOR  NORMAL  OPERATION 
RA( 1 )=ROS/AS*< l.+FMANY*CHANGE) 

DO  170  L=1,LL 

S=(l.-1.574*RA(L)+.35  3*RA(L)**2)/( 1 .-2 . 384*RA ( L ) +1 . 30*RA ( L ) **2 ) 

SS=S-1. 
B=RA(L)*AS/(1.-.5*SS+.486*SS**2-.476*SS**3+.378*SS**4-,158*SS**5) 

POT(L)  =  S*.102276E-23/B**2 
RN(L)  =B 

170  RA(L+1)  =  RA(L)  -  CHANGE  *R0S/AS 
IF(M0DE  -  1)  190,192  ,197 

197  PRINT  196 

196  FORMAT  (13X,7H  ENERGY , 5X , 12H  WAVE  NUMBER , 3X» 20H  WAVE  NUMBER  SQUARE 
ID) 

PRINT  195  ,(E(I),  FK(  I ) ,FK2( I )  ,1=2,11  ) 
19  5  FORMAT  (3E20.4)  , 

PRINT  171 

171  F0RMAT(34H  EFFECTIVE  RANGE/SCATTERING  LENGTH) 

PRINT  172,(RA(L)  L=1,LL) 

172  F0RMAT( E34.10) 
PRINT  194 

194  FORMAT  (13X,7H  RADIUS) 

PRINT  193,  (RN(L),  L=1,LL) 
193  FORMAT  (E20.4) 
PRINT  1192 
1192  FORMAT  ( 16H  POTENTIAL  DEPTH) 

PRINT  2192, (POT( J) ,J=1»JJ) 
2192  FORMAT(F16.2) 
192  PRINT  191 
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191  FORMAK  34H  RADIUS  AND  ENERGY  VALUES  COMPUTED) 
190  DO  187  1=2,11 

ARGB      =  FK(  I  )  *  REND 

CALL  BESSEL<  ARGB , 2 .VALUE , 5 .  0) 

VALJ0( I )  «  VALUE(O) 

VALJK  I  )  =  VALUE(  1  ) 

VALJ2( I )  =  VALUE<2) 

CALL  BESSEL  ( ARGB , 2 .VALUE ,  6.0) 
1  ;  VALNO(  I )  =  VALUE(O) 

VALNK  I  )  =  VALUE(  1  ) 

VALN2( I )  =  VALUE(2) 

IF(MODE  -  1)  187,187,189 
189  PRINT  188,1 
188  FORMAT( 1415) 
187  CONTINUE 

IF(MODE  -  1)  184,186,186 
186  PRINT  185 

185  FORMATC  26H  BESSEL  FUNCTIONS  COMPUTED) 
184  RETURN 

END 
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GAUSSIAN  POTENTIAL 


SUBROUTINE  DERIV  (XCYODY) 

DIMENSION  E(100),FK(100)»FK2(100),RN(100) ♦VALJO(IOO) ,VALJ1( 100) 
1  ,VALJ2(  100)  ,  VALNO(IOO),  VALNl ( 100 ) ,VALN2 ( 100 ) »  RdOOOJt  Y(IOOO), 
2YDOT(1000),  COMPAR0( 100) »  COMPAR 1 ( 100 ) ,  YL ( 1000 ), YLDOT ( 1000 ) »POT ( 1 
3  00)  ,DELTAO( 100)  ,DELT Al ( 100 ) ,XSECO ( 100 )  ,XSEO (100) , TOTAL { 100)  ,RA(10 
40) »FNORM( 100) , BETTER ( 10 » 100 ) 

DIMENSION  YC(30),  DY(30) 

COMMON  E,EINT,EINC,  FKt  FK2,  RINT,  RINC,  RN,  A,  II,  LL,  NN,,  MODE, 
1REND,VALJ0,VALJ1,VALJ2,VALN0,VALN1,VALN2,JJ,M,  NUMBER,  R»I ,L,V,K, 
2N,Y,  YDOT,  YL,  YLDOT,  III,  COMPARO,  C0MPAR1 » POT , VI  NT , V  I NC , MOMENTA , 
3DELTAO,DELTA1,XSECO,XSEC1,TOTAL,RA,ROS,AS,CHANGE,MANY,S,B,FNORM,BE 
ATTER 

DY(1)  =YC(2) 

DY(2)  =  -  (FK2(K)  +  V  *  EXPF(-(XC/A)  **  2))  *  YC(1) 

RETURN 

END 


SUBROUTINE  DER I  VI < XC , YC , DY ) 

DIMENSION  E(100),FK(100),FK2(100),RN(100) ,VALJ0(100)  ,VALJ1 (100) 
1,VALJ2( 100) ,  VALNO(IOO),  VALNl ( 100 ) ,VALN2 ( 100 ) ,  R(IOOO),  Y(1000), 
2YDOT(1000),  COMPAR0(  100) ,  C0MPAR1 ( 100 ) ,  YL ( 1 000 ), YLDOT ( 1000 ), POT ( 1 
300) ,DELTAO( 100)  ,DELTA1 (100) ,XSEC0(100) ,XSEC1 (100) , TOTAL ( 100)  ,RA(10 
40) ,FNORM( 100) , BETTER ( 10,100) 

DIMENSION  YC(30) ♦  DY(30) 

COMMON  E,EINT,EINC,  FK,  FK2,  RINT,  RINC,  RN,  A,  II,  LL ,  NN ,  MODE, 
1REND,VALJ0,VALJ1,VALJ2»VALN0,VALN1,VALN2,JJ,M,  NUMBER,  R,I ,L,V,K, 
2N,Y,  YDOT,  YL,  YLDOT,  III,  COMPARO,  C0MPAR1 , POT ,V I  NT ,V INC , MOMENTA, 
3DELTAO,DELTA1,XSECO,XSEC1,TOTAL,RA,ROS,AS,CHANGE,MANY,S,B,FNORM,BE 
4TTER 

DY( 1  )  =  YC(2) 

DY(2)  =  -(FK2(K)  +  V  *  EXPF ( -< XC/A ) **2 )  -  2. /XC**2 ) *YC ( 1 ) 

RETURN 

END 
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EXPONENTIAL  POTENTIAL 


SUBROUTINE  DERIV  (XOYCDY) 

DIMENSION  E( 100) »FK( 100) »FK2(100) ,RN(100) >VALJ0(100) »VALJ1( 100) 
1 ,VALJ2( 100) ,  VALNO(IOO),  VALN1 ( 100 ) ♦ VALN2 ( 100 ) »  R(1000)t  Y(IOOO). 
2YDOTU000),  COMPAR0(  100)  »  COMP  AR1  (  100  )  ♦  YL ( 1 000 )» YLDOT ( 1000 ) >POT < 1 
300)  ,DELTAO( 100 )  ,DELT Al ( 100 ) ,XSEC0(100) ,XSEC1 (100) »TOTAL(100) >RA(10 
.40) ,FNORM(100)  »  BETTER  (  10  ,  100  ) 

DIMENSION  YC(30) ,  DY(30) 

COMMON  E,EINT,EINC,  FK,  FK2,  RINT,  RINC,  RN,  A.  II,  LL ,  NN,  MODE, 
1REND,VALJ0,VALJ1,VALJ2»VALN0>VALN1,VALN2»JJ»M,  NUMBER*  R»I ,L,V,K, 
2N,Y,  YDOT,  YL,  YLDOT  ♦  I  I  I •  COMPARO.  C0MPAR1 , POT ,V INT , VI NC .MOMENTA , 
3DELTA0 • DELTA1 »XSECO » XSEC1 » TOTAL tRA.ROS. AS »CHANGEt MANY »S»BtFNORMt BE 
4TTER 

DY( 1 )  =  YC(2 ) 

DY(2)  =  -  (FK2(K)  +  V  *  EXPF(-      XC/A ) )  *  YC(1) 

RETURN 

END 

SUBROUTINE  DERI  VI ( XC , YC ,DY ) 

DIMENSION  E( 100) ,FK( 100) ,FK2(100) ,RN(100) ,VALJ0(100) ,VALJ1( 100) 
1,VALJ2( 100) ,  VALNO(IOO),  VALN1 ( 100 ) , VALN2 ( 100 ) ,  R(1000),  Y(1000), 
■2YDOT(1000) ,  COMPAR0( 100) ,  C0MPAR1 ( 100 )  ,  YL ( 1000 ), YLDOT ( 1000 ),  POT ( 1 
300) ,DELTAO( 100 ) ,DELT Al ( 100 ) ,XSECO( 100) ,XSEC1 (100) ,TOTAL(100) ,RA(10 
40)  ,FNORM( 100)  ,BETTER( 10,100) 

DIMENSION  YC(30) ,  DY(30) 

COMMON  E,EINT,EINC,  FK,  FK2,  RINT,  RINC,  RN,  A,  II,  LL ,  NN ,  MODE, 
1REND,VALJ0,VALJ1,VALJ2»VALN0»VALN1»VALN2,JJ»M,  NUMBER,  R,I,L,V,K, 
2N,Y,  YDOT,  YL,  YLDOT,  III*  COMPARO,  C0MPAR1 , POT , VI  NT , VI NC, MOMENTA, 
3DELTAO,DELTA1,XSECO,XSEC1,TOTAL,RA,ROS,AS,CHANGE,MANY,S,B,FNORM,BE 
4TTER 

DY( 1)  =  YC(2) 

DY(2)  =  -  (FK2(K)  +  V  *  EXPF(-      XC/A)  -2,  /XC**2 ) *YC ( 1 ) 

RETURN 

END 
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YUKAWA  POTENTIAL 


SUBROUTINE  DERIV  (XOYCDY) 

DIMENSION  E(100),FK(100),FK2(100),RN(100)  »VALJ0(100)  ,VALJ1(  100) 
1 »VALJ2( 100)  ,  VALNO(IOO),  VALN1 (  100 )  * VALN2 ( 100  )  ♦  R(1000)»  Y(1000)t 
2YDOTCL000)  ,  COMPARO ( 1 00 )  »  COMPARl ( 100 )  ,  YL ( 1 000 ), YLDOT ( 1000 ), POT ( 1 
3  00)  ,DELTAO( 100)  ,DELT Al ( 100  )  ,XSECO ( 100 )  ,XSEC1 (100)  , TOTAL (100)  ,RA(10 
40)  ,F:NORM(  100)  ,  BETTER  (  1  0  ,  100  ) 

DIMENSION  YC(30) ,  DY(30) 

COMMON  E,EINT,EINC,  FK,  FK2,  RINT,  RINC,  RN,  A,  II,  LL,  NN ,  MODE, 
1REND,VALJ0,VALJ1,VALJ2,VALN0,VALN1»VALN2»JJ»M,  NUMBER,  R,I ,L,V,K, 
2N,Y,  YDOT,  YL,  YLDOT  ,  III,  COMPARO,  C0MPAR1 , POT , VI  NT , VI NC , MOMENTA , 
3DELTAO,DELTA1,XSECO,XSEC1,TOTAL,RA,ROS,AS,CHANGE,MANY,S,B,FNORM,BE 
ATTER  |  . 

DY( 1 )  =  YC(2)  4 

DY(2)  =  -  (FK2CK)  +  V  *  ( EXPF ( -XC/A ) ) / ( XC/A ) ) *YC ( 1 ) 

RETURN 

END 

SUBROUTINE  DER I  VI ( XC , YC ,DY ) 

DIMENSION  E(100),FK(100) »FK2( 100 ) , RN ( 100 ) , VAL JO ( 100 )  , VALJ1 ( 100 ) 
1 ,VALJ2( 100) ,  VALNO(IOO),  VALN1 ( 100 ) , VALN2 ( 100 ) ,  R(IOOO),  Y(IOOO), 
2YDOT(1000),  COMPAR0( 100) ,  COMPARl ( 100 ) ,  YL ( 1 000 ), YLDOT ( 1 000 ), POT ( 1 
3  00)  ,DELTAO( 100 ) ,DELT Al ( 100 ) ,XSEC0(100) ,XSEC1 (100) , TOTAL ( 100)  ,RA(10 
40) ,FNORM( 10  0) ,BETTER( 10,100) 

DIMENSION  YC(30),  DY(30) 

COMMON  E,EINT,EINC,  FK ,  FK2,  RINT,  RINC,  RN,  A,  II,  LL ,  NN ,  MODE, 
1REND,VALJ0,VALJ1,VALJ2,VALN0»VALN1 ,VALN2,JJ»M,  NUMBER,  R,I ,L,V,K, 
'2N,Y,  YDOT,  YL,  YLDOT,  III,  COMPARO,  COMPARl , POT ,V INT  ,VINC , MOMENTA , 
3DELTAO,DELTA1,XSECO,XSEC1»TOTAL,RA,ROS,AS,CHANGE,MANY»S,B,FNORM»BE 
4TTER 

DY( 1 )  =  YC(2) 

DY(2)  =  -  (FK2(K)  +  V  *  (  EXPF  ( -XC/A ))/( XC/A )  -2./XC'**2)  *  YC(1) 

RETURN 

END 
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SQUARE  WELL  POTENTIAL 


SUBROUTINE  DERIV  (XCYODY) 

DIMENSION  E( 100) ,FK( 100) »FK2( 100) ,RN(100) ,VALJ0(100) »VALJ1( 100) 
1 ,VALJ2( 100)  ,  VALNO(IOO)*  VALN1 ( 100 ) , VALN2 ( 100 ) »  R(IOOO).  Y(IOOO). 
2YD0T(1000),  COMPAR0( 100) .  COMPARl ( 100 ) .  YL ( 1 000 ) »YLOOT ( 1000 ) »POT ( 1 
3  00)  ,DELTAO( 100)  ,DELTA1 (100)  ,XSECO( 100)  ,XSEC1 (100) » TOTAL (100)  »RA(10 
40) ,FNORM( 100) , BETTER ( 10 , 100 ) 

DIMENSION  YC(30) ,  DY(30) 

COMMON  E,EINT,EINC,  FK  ,  FK2,  RINT,  RINC,  RN,  A»  II,  LL,  NN ,  MODE, 
1REND,VALJ0,VALJ1,VALJ2»VALN0,VALN1 ,VALN2»JJ»M,  NUMBER,  R» I »L,V,K, 
2N,Y,  YDOT,  YL,  YLDOT ,  III,  COMPARO,  COMPARl ,POT ,V I  NT , VI NC , MOMENTA, 
3DELTAC,DELTA1,XSECO,XSEC1,TOTAL,RA,ROS,AS,CHANGE,MANY,S,B,FNORM,BE 
4TTER 

DY(1)  =YC(2) 

IF(XC  -A)  1,1,2 
.  DY(2  )  =  -(FK2(K)  +  V)  *  YC ( 1 ) 

RETURN 
:  DY(2)  =-<FK2(K)  )*  YC ( 1  ) 

RETURN 

END 


SUBROUTINE  DERI  VI ( XC , YC  ,DY  ) 

DIMENSION  E( 100) ,FK( 100) ,FK2( 100) ,RN( 100) ,VALJ0(100) ,VALJ1( 100) 
1 ,VALJ2( 100) ,  VALNO(IOO),  VALN1 ( 100 ) , VALN2 ( 100 ) ,  R(IOOO),  Y(1000), 
2YDOT(1000),  COMPAR0( 100) ,  COMPARl ( 100 ) ,  YL ( 1000 ), YLDOT ( 1000 ), POT ( 1 
3  00) ,DELTAO( 100 ) ,DELT Al ( 100 ) ,XSEC0(100) ,XSEC1 (100) , TOTAL ( 100) ,RA(10 
40) ,FNORM( 100) ,BETTER( 10,100) 

DIMENSION  YC(30) ,  DY(30) 

COMMON  E,EINT,EINC,  FK,  FK2,  RINT,  RINC,  RN,  A,  II,  LL ,  NN ,  MODE, 
1REND,VALJ0,VALJ1,VALJ2,VALN0,VALN1,VALN2,JJ»M,  NUMBER,  R,I ,L,V,K, 
2N,Y,  YDOT,  YL,  YLDOT,  III,  COMPARO,  COMPARl , POT ,V I  NT ,V I NC , MOMENTA , 
3DELTAO,DELTA1,XSECO,XSEC1,TOTAL,RA,ROS,AS,CHANGE,MANY,S,B,FNORM,BE 
4TTER 

DY(1)  =YC(2) 

IF (XC  -A)  1,1,2 
.  DY(2)  =  -(FK2(K)  +  V  -  2 ./XC**2 )*YC ( 1 ) 

RETURN 
!  DY(2)  =  -(FK2(K)  -  2./  XC**2 )  *YC(1) 

RETURN 

END 
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THE  FOLLOWING  DEFINITIONS  FOR  B  MUST  BE  SUBSTITUTED  IN  READY2  FOR  VARYING 


CATTERIN6  LENGTH.  THE  DEFINITION  PROVIDED  IN  READY2  IN  THE  LISTING  IS  FOR 


ARIATIONS  IN  EFFECTIVE  RANGE. 


YUKAWA  WELL 


B=ROS/(1.-1.369*SS+1.0  9  3*SS**2-1.127*SS**3+1.00  5*SS**4-.324*SS**5) 


SQUARE  WELL 


B=ROS/(l.-.5*(S-l. ) +.486* (S-l. ) **2-.476* ( S-l . ) **3+. 3  7  8* ( S-l . ) **4 


~.158*(S-1. )**5 ) 


EXPONENTIAL  WELL 


B=ROS/( l.-.904*SS+.74  5*SS**2-.543*SS**3+.174*SS**4) 


GAUSSIAN  WELL 


B=ROS/( l.-.641*SS+.568*SS**2-.721*SS**3+1.172*SS**4-1.015*SS**5) 
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APPEND  !X  IV 
SAMPLE  RESULTS 
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These  results  are  for  a  yukawa  potential  with  r  *  f\$ed   at 

o,s 

2.46  Fm  and  a  varying  In  I  16%  deviation  from  a  *  =  23.677  FM. 
s  s 
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RACTUS      POTENTIAL  CEPTH 
.1Q86E-12  .5621E+02 

ENERGY  PHASE  ANGLE  CRCSS  SECTION, 
L=0         L=0 


.10 

.20 

.30 

.40 

.50 

.60 

.70 

.80 

.90 

1.00 

1.10 

1.20 

1.30 

1  .40 

1.50 

1.60 

1.70 

1  .80 

1.90 

2.00 

2.10 

2.20 

2.30 

2.40 

2.50 

2.60 

ENERGY 

.10 

.20 

.30 

.40 

.50 

.60 

.70 

80 

90 

00 

10 

20 

30 

40 

50 

60 

70 

80 


1.90 
2.00 
2.10 
2.20 
2.30 
2.40 
2.50 
2.60. 


9113 
1037 
1094 
1125 
1145 
1158 
1166 
1172 
1176 
1178 
1180 
1180 
1180 
1180 
1  179 
1178 
1176 
1175 
1173 
1172 
1  170 
1168 
1166 
1164 
116  1 
1159 
.NO 


E  +  00 

E  +  01  . 

E+01 

E  +  01 

E+0  1 

E  +  01 

E+01 

E  +  01 

£  +  01 

E+01 

E  +  01 

E  +  01 

E  +  01 

E+01 

E+01 

t+01 

E  +  01 

E  +  01 

E+01 

E  +  01 

E+01 

E  +  01  - 

E+01 

E  +  01 

E  +  01 

E+01 

REALIZE 

1  136696 

1036931 

1065E04 

1C54102 

10*46656 

1041 494 

10377CC 

1034789 

1032432 

10  3060c 

1029C48 

102773? 

1026607 

1025630 

1024773 

1024016 

1023340 

1022733 

1022185 

1021686 

1021231 

1020812 

1020427 

1020C7C 

1019738 

10JL9429 


.3257E 
.1933E 
.1372E 
.1062E 
.8649E 
.7290E 
.6295E 
.5535E 
.4936E 
.4451E 
.4051 E 
.3716E 
.3430E 
.3184E 
.2970E 
.278 1  £ 
.26  15E 
.2466E 
.2333E 
.2213E 
.21C5E 
.2006E 
915E 
832E 
755E 
685E 
CROSS  S 


C 

7E  +  C 
1E+C 
7E  +  C 
OE+0 
4E  +  C 
7E  +  C 
2E  +  0 
2L  +  0 
1E  +  C 
2E  +  C 
9E  +  0 
8E  +  0 
1E  +  C 
CE  +  C 
6E  +  0 
1E  +  0 
5E  +  C 
8E  +  C 
3E+0 
8E  +  0 
1E  +  C 
9E  +  0 
2E+C 
2E+0 
6E  +  C 
6E  +  0 


-22 

-22 

-22 

-22 

-23 

-23 

-23 

-23 

-23 

-23 

-23 

-23: 

-23 

-23' 

-23   ■ 

-23 

-23 

-23 

-23 

-23 

-23 

-23 

-23 

-23 

-23 

-23 

ECTION 


92 


EFFECT  I  VTT"R AWE /"SCATTERTKXT  LENGTH 

-.9143050217E-01 
RADIUS      POTENTIAL  DEPTH 
.1C83E-12  .56141 E+02 

ENERGY  PHASE  ANGLE  CRGSS  SECTION,' 
L=0         L=0 


.10 

.20 

.30 

.UO 

.50 

.60 

.70 

.80 

.90 

1.00 

1.10 

1.20 

1.30 

1.40 

1.50 

1.60 

1.70 

1  .80 

1.90 

2.00 

2.10 

2.20 

2.30 

2.40 

2.50 

2.60 

ENERGY 

.10 

.20 

.30 

.40 

.50 

.60 

.70 

.80 

.90 

1.00 

1.10 

1.20 

1.30 

1.40 

1.50 

1.60 

1.70 

1.80 

1.90 

2.00 

'2.10 

2.20 

2.30 

2.40 

2.50 

2.6CL 


6913E 

020E 

079E 

1  12E 

133E 

146E 

156E 

162E 

166E 

169E 

171E 

172E 

172E 

172E 

172E 

171E 

170E 

168E 

167E 

1o5E 

164E 

162E 

160E 

158E 

156E 

154E 

NOR 

.  1 

.1 

.1 

.1 

.1 

.1 

.1 

.1 

.1 

.1 

.1 

.1 

.1 

.1 

.1 

.1 

.1 

.1 

.1 

.1 

.1 

.1 

.1 

.1 

.1 

__  .1 


+  00 
+  0 
+0 
+  0 
+  0 
+  0 
+  0 
+  0 
+  0 
+0 
+  0 
+  0 
+  0 
+  0 
+  0 
+  0 
+  0 
+  0 
+  0 
+  0 
+  0 
+0 
+  0 

■5-0 

+0 

+0 

iVALIZ 

101  14 

06497 

04936 

04065 

03509 

03122 

02837 

02618 

024U4 

02303 

02185 

02086 

02001 

01927 

01863 

01805 

01754 

01708 

01667 

01629 

01595 

01563 

01534 

01507 

01482 

OlU±5S 


•  o 


ED  C 

34E  + 
39E  + 
S6E  + 
9  1E  + 
36E  + 
U3E  + 
43E  + 
147E  + 
73E  + 
35E  + 
9CE  + 
66E  + 
61E  + 
82E  + 
13E  + 
89E  + 
84E  + 
98E  + 
53E+ 
814E  + 
39E  + 
77E  + 
6  1E+ 
62E  + 
55E  + 
.LS£  + 


155 
.1894 
.1351 
.1048 
.8553 
.7218 
.6238 
.5489 
.4897 
.4419 
.4023 
.3691 
.3403 
.316*4 
.2952 
.2765 
.2600 
.2453 
.2321 
.2202 
.2094 
.1995 
.1906 
.1823 
.1747 
.1677 
ROSS 
01 
01 
01 
01 
01 
01 
01 
01 
01 
01 
01 
01 
01 
01 
01 
01 
01 
01 
01 
01 
01 
01 
01 
01 
01 
CI 


E-22 
E-22 
E-22 

E-22 

E-23 

E-23 

E-23 

E-23 

E-23 

E-23 

E-23 

E-23 

E-23 

E-23 

E-23 

E-23 

E-23 

E-23 

E-23 

E-23 

E-23 

E-23 

E-23 

E-23 

E-23 

E-23. 

SECTICN 


33 


EFFECTIVE  RANGE /SCTnTTERTNG  LENGTH" 

-♦9558643409 E-01 
RADIUS      POTENTIAL  DEPTH 
.1C8GE-12  .5661E+C?  . 

ENERGY  PHASE  ANGLE  CRCSS  SECTICN, 
L  =  C         L=0 


.10 

.20 

.30 

.1+0 

.50 

.60 

.70 

.80 

90 

00 

10 

20 

30 

1  .40 

1.50 

1.60 

1.70 

1  .80 

1.90 

2.00 

2.10 

.20 

30 


40 
50 
60 


ENERGY 

.  10 

.20 

.30 

.40 

.50 

.60 

.70 

.80 

.90 

1.00 

1.10 

1.20 

1.30 

1.40 

1.50 

1.60 

1.70 

1.80 

1  .90 

2.00 

2.10 

2.20 

2.30 

2.40 

2.50 

2.60 


8719 
1004 
1064 
1099 
1121 
1135 
1145 
1152 
1157 
1161 
1163 
llo4 
1165 
1165 
1165 
1164 
1163 
1162 
1161 
1159 
1158 
1156 
1154 
1152 
1150 
1148 
NO 


E  +  00 

E+0  1 

E  +  01 

E+01 

E  +  01 

E+0  1 

£+01 

E+01 

E  +  01 

E  +  01 

E  +  01 

E  +  01 

E  +  01 

E  +  01 

E  +  01 

E+01 

E+Cl 

E  +  01 

E+01 

E+01 

E+01 

E  +  01 

E  +  01 

E  +  01 

E+01 

E+01 

RiVAL 

1066 

1043 

1032 

1027 

1023 

1020 

1018 

1017 

1016 

1015 

1014 

1013 

1013 

1012 

1012 

1012 

1011 

1011 

1011 

1010 

1010 

1010 

1010 

1010 

1009 

1009 


IZEC 
4810 
138C 
903  C 
1472 
452? 
8768 
9756 
5126 
3506 
4041 
6173 
9521 
3817 
8867 
4526 
0684 
7257 
4178 
1394 
8863 
6549 
4425 
2467 
065U 
897C 
74CC 


.3056E 
.1855E 
.1329E 
.1034E 
.8457E 
.7145E 
.6181 E 
.5443E 
.4859E 
.4386E 
.3994E 
.3666E 
.3386E 
.3144E 
.2934E 
.2749E 
.2585E 
.2439E 
.2  308E 
.2190E 
.2083E 
985E 
896E 
814E 
733E 
669E 
CROSS  S 

E  +  0 

E  +  C 

E  +  C 

E  +  C 

E  +  C 

E  +  C 

E+0 

E  +  0 

E  +  C 

E  +  C 

E+0 

E+0 

E  +  0 

E  +  C 

E+0 

E+0 

E  +  C 

E  +  C 

E+0 

E+0 

E  +  C 

E  +  0 

E  +  0 

E  +  0 

E  +  C 

E  +  C 


-22 

-22 

-22 

-22 

-23 

-23 

-23 

-23 

-23 

-23 

-23 

-23 

-23 

-23 

-23 

-23 

-23 

-23 

-23 

-23 

-23 

-23 

-23 

-23 

-23 

-23  -—■ 

ECTICN 


9<* 


EFFECTIVE  RANGE/SCATTERING  LENGTH 

-.9974236601E-01 
RADIUS      POTENTIAL  DEPTH 
.1C77E-12  .5680E+C2 

ENERGY  PHASE  ANGLE  CRCSS  SECTICN, 
L=0         L=C 


.10 
.20 
.30 

.140 

.50 

.60 

.70 

.80 

.90 

1.00 

1.10 

1.20 

1.30 

1.40 

1.50 

1.60 

1.70 

1  .80 

1.90 

2.00 

2.10 

2.20 

2.30 

2.40 

2.50 

2.60 

ENERGY 

.10 

.20 

.30 

.40 

.50 

.60 

.70 

.80 

.90 

1.00 

1.10 

1.20 

1.30 

1.40 

1.50 

1.60 

1.70 

1.80 

1.90 

2.00 

2.10 

2.20 

2.30 

2.40 

2.50 

2.60. 


8531 
9875 
1050 
1086 
11C9 
1125 
1135 
1143 
11U8 
1152 
1155 
1156 
1157 
1158 
1157 
1157 
1156 
1155 
1154 
11-3 
1152 
1150 
1148 
11^7 
1145 
111)3 
NO 


E  +  00 

E+00 

E+01 

E  +  01 

E+0  1 

E+01 

E  +  01 

E  +  01 

E  +  01 

E+01 

E  +  01 

E  +  01 

E  +  01 

E+01 

E+01 

E+01 

E+01 

E+01 

E  +  01 

E  +  01 

E+01 

E  +  01 

E+01 

E  +  0  1 

E+01 

E+0  1 

RMALIZE 

1032755 

1021467 

1016437 

1013587 

10117149 

1010464 

1009514 

100878  1 

1008199 

1007724 

1007329 

1006995 

1006708 

1006459 

100624  1 

1006047 

1005875 

100572C 

100558C 

1005452 

1005336 

1005229 

1005130 

1005059 

1004954 

1004875 


.2959E 
. 18  17  E 

.1308E 
.1021E 
.8360E 
.7073E 
.6124E 
.5396E 
.4820E 
.4352E 
.3966E 
.3641 E 
.3363E 
.3124E 
.2916E 
.2733E 
.2570E 
.2425E 
.2296E 
.2178E 
.2072E 
975E 
886E 
8C5E 
730E 
661E 
CROSS  S 


D 

5E  +  0 
1E  +  0 
6E+C 
4E  +  0 
8E  +  0 
7E  +  C 
2E+0 
7E  +  0 
2E  +  0 
3E  +  0 
2E  +  0 
OE  +  0 
3E  +  0 
4E  +  C 
1E  +  0 
8E  +  0 
3E  +  C 
3E  +  C 
2E+0 
7E  +  0 
3E  +  C 
5E  +  C 
7E  +  0 
4E  +  0 
6E  +  0 
6E  +  G 


-22 

-22 

-22 

-22 

-23 

-23 

-23 

-23 

-23 

-23 

-23 

-23 

-23 

-23 

-23 

-23 

-23 

-23 

-23 

-23 

-23 

-23 

-23 

-23  . 

-23 

-23-  . 

ECTION 


95 


EFFECTIVE  RANGE/SC  ATTER11TG~L~ET\~GTF 

-.1C38982979E+00 
RADIUS      POTENTIAL  CEPTF 
.1C74E-12  .5700E+C2 

ENERGY  PHASE  ANGLE  CRCSS  SECTICN, 
L=0         L=0 


.10 

.20 

.30 

.40 

.50 

.60 

.70 

.80 

.90 

1  .00 

1.10 

1.20 

1.30 

1.40 

1.50 

1.60 

1.70 

1.80 

1  .90 

2.00 

2.10 

2.20 

2.30 

2.40 

2.50 

2.60 

ENERGY 

.10 

.20 

.30 

.40 

.50 

.60 

.70 

.80 

90 

00 

10 

20 

30 

1.40 

1.50 

1.60 

1.70 

1.80 

1  .90 

2.00 

2.10 

2.20 

2.30 

2.40 

2.50 

2.60 


8350E+00 
9717E+00 
1036E+G 
074L+0 
C98E+0 
1  14E  +  0 
125E+0 
153E+0 
139E+0 
143E+0 
146E+0 
148E+0 
150E+0 
150E+0 
150E+0 
150E+0 
150E+0 
149E+0 
148E+0 
147E+0 
146 E+0 
144E+0 
143E+0 
14  1 E+0 
140 E+0 
138 E+0 
NORMAL 
.  1000 
.  1000 
.1000 
.1C00 
.  1000 
.  1000 
.1000 
•  1000 
.1000 
.1000 
.1000 
.1000 
.  1C00 
.  1000 
.1000 
.1000 
.1000 
.1000 
.1000 
.1000 
.  1000 
.  1000 
.1000 
.  1000 
.1000 
.1000 


.2 

.1 
.1 
.1 

.8 
.6 
.6 
.5 
.4 
.4 
.3 

.  -' 

.3 
.3 
.2 

.2 
.2 
.2 
.2 
.2 


IZEC  CRD 

OOOOE+0 

OOCCE+0 

OOOCE+C 

OOOOE+0 

COOOE+O 

OOOCE+C 

COCCE+0 

OOOOE+O 

OOOOE+0 

OOCCE+C 

COCCE+0 

OOOOE+0 

OOOOE+C 

OOCCE+C 

CCCCE+O 

OOOOE+0 

OOOOE+0 

OOCCE+0 

COCCE+0 

OOOOE+0 

OOCCE+C 

OOCCE+0 

CCCCE+O 

OOCOE+0 

OOCCE+0 

OOC_CE+.0 


866E 
778E 
287E 
007E 
263E 
999E 
066E 
349E 
781E 
319E 
937E 
615  E 
341  E 
104E 
898E 
716E 
555E 
412E 
283E 
166E 
061E 
965E 
877  E 
796E 
721  E 
653E 

ss  s 


-22 

-22 

-22   • 

-22 

-23 

-23 

-23 

-23 

-23 

-23 

-23 

-23 

-23 

-23 

-23 

-23 

-23 

-23 

-23 

-23 

-23 

-2  3 

-23 

-23 

-23 

-2*  - 

ECT  Toisf 


96 


EFFECTIVE  RANGE/SCATTERING  LENGTH 

-.1C8C542298E+00 
RADIUS      POTENTIAL  DEPTH 
.1C71E-12  .5719E+C2 

ENERGY  PHASE  ANGLE  CRCSS  SECTION, 
L=0   ■      L=0 


.10 
.20 
.30 

.uo 

.50 

.60 

.70 

.80 

.90 

1.00 

1.10 

1.20 

1.30 

1  .40 

1.50 

1.60 

1.70 

1  .80 

1.90 

2.00 

2.10 

2.20 

2.30 

2.40 

2.50 

2.60 

ENERGY 

.10 

.20 

.30 

.40 

.50 

.60 

.70 

.80 

.90 

1.00 

1.10 

1.20 

1.30 

1.40 

1.50 

1.60 

1.70 

1  .80 

1.90 

2.00 

2.10 

2.20 

2.30 

2.U0 

2.50 

2.60 


8175 
9562 
1022 
1061 
1086 
1103 
1115 
1124 
1  130 
1135 
1138 
1141 
1142 
1143 
1144 
1144 
1143 
1143 
1  142 
114  1 
1140 
1  139 
1137 
1136 
1134 
1133 
NO 


E  +  00 

E  +  00 

E+0 

E  +  0 

E  +  0 

E+0 

E  +  0 

E  +  0 

E+0 

E+0 

E+0 

E  +  0 

E  +  0 

E  +  0 

E  +  0 

E+C 

E  +  0 

E+0 

E  +  0 

£  +  0 

E  +  0 

E  +  0 

E  +  0 

E  +  0 

E+0 

E  +  0 

RHALIZ 

968236 

978771 

983616 

986403 

988218 

989494 

990442 

991 175 

991759 

992237 

992634 

992971 

993260 

993512 

993732 

993928 

994102 

994259 

9944C1 

994530 

994648 

994756 

994856 

994948 

995034 

9951  14 


.2775E 
.1741E 
.1266E 
.9934E 
.8166E 
.6926E 
.6008E 
.53C2E 
.4741E 
.4286E 
.3908E 
.3590E 
.3318E 
.3084E 
.2880E 
.2700E 
.2540E 
.2398E 
.2270E 
.2  155  E 
.2050E 
.195ME 
.1867E 
.1787E 
.1713E 
.1645E 

ED  CROSS  S 

62E+C0 

50E+C0 

C2E+C0 

98E+00 

1  1E  +  C0 

47E+C0 

57E+00 

5  4 E+00 

86E+GO 

16E+C0 

82E+00 

62E+00 

84E+C0 

16E+C0 

77E+00 

18E+00 

61E+C0 

41E+00 

25E+C0 

26E+00 

2CE+C0 

52E+C0 

4  1  E+00 

88E+00 

79E+C0 

85E+00 


-22 

-22 

-22 

-23 

-23 

-23 

-23 

-22 

-23 

-23 

-23 

-23 

-23 

-23' 

-23 

-23 

-23 

-23 

-23 

-23 

-23 

-23 

-23 

-23 

-23 

-23 

ECTICN 


97 


EFFECTIVE  RANGE/SCATTERING  LENGTH 

-.1  1221C1618E+0G 
■RACIUS      POTENTIAL  CEPTh 
.1C68E-12  .5738El+£2 

ENERGY  PHASE  ANGLE  CRCSS  SECTION, 
.  L  =  0   .      L=0  I 


.10 

.80C5E+00 

•  2686E- 

-22 

.20 

.9U12E+C0 

.1703E- 

-22 

.30 

.  1009E+0  1 

.121+5E- 

-22 

'  .UO 

.  1049E+01 

.9797E- 

-23 

.50 

.  1075E+01 

.8C68E- 

-23 

.60 

.  1093E  +  01 

.6852E- 

-23 

.70 

.  1  1C5E+C1 

•  5950E- 

-23 

.80 

.1  1  15E+01 

.525UE- 

-23 

.'90 

.  1  122E+0  1 

.U7C2E- 

-23 

1  .00 

.  1127E+01 

.1+2525- 

-23 

1.10 

,  1  130E+C1 

.3879E- 

-23 

, 

1.20 

,  1  133E+0  1 

.3561+E- 

-23 

1.30 

.  1135E+01 

.3296E- 

-23 

1  .1+0 

.  1136E+01 

.3061+E- 

-23 

1  .50 

. 11 37E+01 

.2861E- 

-23 

1.60 

, 1137E+01 

.2683E- 

-23 

1.70 

. 1137E+C1 

.2525E- 

-23 

1.80 

.11360+01 

.2381+E- 

-23 

1.90 

1136E+01 

.2257E- 

-23 

2.00 

, 1 135E+01 

.2143E- 

-23 

2.10 

. 1 15LE+01 

.2029E- 

-23 

2.20 

. 11 33E+01 

.  19UUE- 

-23 

2.30 

1152E+01 

.1857E- 

-23 

2.1+0 

. 1 130E+C1 

.1778E- 

-23 

2.50 

, 1129E+01 

.1704E- 

-23 

2.60 

, 1127E+01 

NORMALIZED 

.  1636E- 
CROSS    Si 

-23 

ENERGY 

ECTICN 

.10 

.937U7776E+C0 

.20 

.95781212E 

+  00 

.30 

.96730933E+G0 

1 

.40 

.9728  1728E 

+  00 

! 
i 

.50 

.9761+1799E  +  00 

.60 

.97895922E 

+  C0 

i 

.70 

.980851C3E+00 

• 

.80 

.98231599E 

+  00 

.90 

.98348532E 

+  C0 

1.00 

. 981+1+1+ 11+3E  +  C0 

1.10 

.98523866E 

+  00 

1.20 

.98591430E 

+00 

1.30 

-   .986U91+80E 

+  00 

1  .40 

. 986999U  i+E  +  CO 

1.50 

.9874426CE+0C 

1.60 

.98783524E 

+  00 

1.70 

.98818582E+00 

1.80 

.988501C1+E 

+  00 

1.90 

.98878622E+00 

; 

2.00 

.98901+565E+00 

1 

2.10 

.9892S285E+00 

! 

2.20 

.9895C072E 

+  00 

2.30 

.9897C165E 

+  00 

2.40 

.98988768E 

+00 

2.50 

.9900605CE 

+  C0 

i 

2.60 

.99022157E+C0 

1 

A 
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EFFECTIVE    RANGE/SCATTERING    LENGTH 

-.1 163660937E+00 
RACIUS  POTENTIAL    DEPTH 

. 1C65E-12  ,5757.SAfi5 

ENERGY    PHASE    ANGLE    CRCSS    SECTION, 
L=0  L=0 


.10 

.7842E+00 

.2601E-22 

.20 

.9264E+00 

.1667E-22 

.30 

.9957E+00 

.1224E-22 

.40 

.1037E+01 

.9661E-23 

.50 

.1064E+01 

.7971E-23 

.60 

.  108  2E+01 

.6778E-23 

.70 

,  1096E  +  01 

.5892E-23 

.80 

.1 1C6E+01 

.5207E-23 

.90 

.11 13E+01 

.4662E-23 

1.00 

.  1118E  +  01 

.4218E-23 

1.10 

,  1122E+01 

.3850E-23 

1.20 

.1125E+01 

.3539E-23 

1.30 

,1 127E+01 

.3273E-23 

1.40 

.1129E+01   " 

.30144E-23 

1.50 

1150E+01 

.2843E-23 

1.60 

.  1 130E  +  01 

.2667E-23 

1.70 

,  1130E+01 

.2510E-23 

1.80 

,  1130E+01 

.2370E-23 

1.90 

, 1130E+01 

.2244E-23 

2.00 

,  1129E  +  01 

.2131E-23 

2.10 

.  1128E+01 

.2028C-23 

2.20 

.1127E+01 

.1934E-23 

2.30 

,  1126E+01 

.1348E-23 

2.40 

.1125E+01 

.1769E-23   ' 

2.50 

.  1124E  +  01 

.1696E-23 

2.60 

,  1122E+01 

.1628E-23 

"VVWV 

ENERGY 

NORMALIZED 

CROSS  SECTION 

( 

•  .10 

.9077276CE+00 

.20 

.9371U849E 

+  00 

.30 

.951101 cSE+CO 

.40 

.95925670E 

+00 

.50 

.96461287E 

+  C0 

.'60 

.96840489E 

+  C0 

. 

.70 

.97123418E 

+  00 

.80 

.97342875E+00 

.90 

.97518273E+00 

1.00 

.97661832E 

+  00 

1 

1.10 

.9778163HE+C0 

1.20 

.97883222E 

+  00 

1.30 

.97970572E 

+  C0 

1.40 

.98046532E 

+  00 

1.50 

.98113263E 

+  00 

1.60 

.98172U05E 

+00 

1.70 

.98225227E+C0 

1.80 

.98272732E+00 

1.90 

.98315716E+00 

2.00 

.93354827E 

+00 

2.10 

.98390590E+CO 

2. 20 

.98423442E 

+  00 

;' 

2.30 

.98453743E+00 

1 

2. 40 

.93481798E+00 

1 

2.50 

.98507864E 

+  00 

2.60 

.9853215SF_+00 

.  j 
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EFFECTIVE  RANGE/SCATTERING  LENGTH. 

-.12C5220256E+00 
RADIUS      POTENTIAL  DEPTH 
.1C62E-12  .5776E+C2 

ENERGY  PHASE  ANGLE  CRCSS  SECTION, 
L=0         L=0 


.10 

.20 

.30 

.40 

.50 

.60 

.70 

.80 

.90 

1.00 

1.10 

1.20 

1.30 

1.40 

1.50 

1.60 

1.70 

1  .80 

1  .90 

2.00 

2.10 

2.20 

2.30 

2. HO 

2.50 

2.60 

ENERGY 

.  10 

.20 

.30 

.40 

.50 

.60 

.70 

.80 

.90 

1.00 

1  .10 

1.20 

1.30 

1.40 

1.50 

1.60 

1.70 

1.80 

1  .90 

2.00 

2.10 

2.20 

2.30 

2.40 

2.50 

2.60 


7684E+00 

9120E+00 

9828E+00 

1025E+0 

1053E+0 

1072E+0 

1086E+0 

1097E+0 

1  1C4E+0 

I  1  10E+0 

II  14E  +  G 
11 18E+0 
1120E+0 
1122E+G 
1 123E+0 
1124E+0 
1124E+0 
1124E+0 
1124E+0 
1123E+0 
1122E+0 
1 122E+0 
1121E+0 
1 1 19E  +  0 
1 1 18E+0 
1117E+0 

NORMAL 
.8789 
.9168 
.9350 
.9457 
.9528 
.9578 
.9616 
.9645 
.9668 
.9637 
.9703 
.9717 
.9728 
.9759 
.9748 
.9755 
.9763 
.9769 
.9775 
.9730 
.9785 
.9789 
.9793 
.9797 
.9800 
.9304 


IZE 
833 
036 
13C 
379 
152 
416 
004 
209 
581 
731 
725 
298 
973 
131 
C59 
97  U 
045 
406 
162 
401 
192 
593 
654 
Ml2 
906 
162 


D 

IE 
3E 
CE 

5E 

c  c 

5E 
5E 
6E 
9E 
7E 
6E 

2E 
6E 
2E 
IE 
3E 
CE 
7E 
4E 
hE 
9E 
IE 
6E 
8E 
7F 


.2519E-22 
.1630E-22 
.1203E-22 
•9525E-23 
.7873E-23 
.67C4E-23 
.5833E-23 
.5159E-23 
.4622E-23 
.41S4E-23 
.3320E-23 
.3513E-23 
.325GE-23 
.3023E-23 
.2825E-23 
.2650E-23 
.2495E-23 
.2356E-23 
.2231E-23 

1  19E-23  i 

017E-23 

923E-23 

838E-23 

760E-23 

637E-23 

620E-23  - 

SS  SECTION 


.2 
.1 
.1 

.1 
.1 
.1 

CRO 
+  00 
+  C0 
+  00 
+  00 
+  00 
+  00 
+  00 
+  00 
+  00 
+  00 
+  00 
+  00 
+  00 
+  00 
+  00 
+  00 
+  00 
+  00 
+  00 
+  00 
+  00 
+  00 
+  00 
+00 
+  00 
+  C0 


*  v^'v^r 
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